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Abstract: We consider the possibility that the total dark energy (DE) of the Universe is 
made out of two dynamical components of different nature: a variable cosmological term, 



in 



A, and a dynamical "cosmon", X, possibly interacting with A but not with matter - which 
remains conserved. We call this scenario the AXCDM model. One possibility for X would 
be a scalar field x, but it is not the only one. The overall equation of state (EOS) of the 
AXCDM model can effectively appear as quintessence or phantom energy depending on the 
mixture of the two components. Both the dynamics of A and of X could be linked to high 
energy effects near the Planck scale. In the case of A it may be related to the running of this 
parameter under quantum effects, whereas X might be identified with some fundamental 
field (say, a dilaton) left over as a low-energy "relic" by e.g. string theory. We find that the 
S^ dynamics of the AXCDM model can trigger a future stopping of the Universe expansion and 

can keep the ratio pd/ Pm (DE density to matter-radiation density) bounded and of order 
1. Therefore, the model could explain the so-called "cosmological coincidence problem". 
This is in part related to the possibility that the present value of the cosmological term can 
be Ao < in this framework (the current total DE density nevertheless being positive). 
However, a cosmic halt could occur even if Ao > because of the peculiar behavior of X 
as "Phantom Matter". We describe various cosmological scenarios made possible by the 
composite and dynamical nature of AXCDM, and discuss in detail their impact on the 
cosmological coincidence problem. 
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1. Introduction 

Modern experimental cosmology offers us a strikingly accurate picture of our Universe 
which has been dubbed the cosmological concordance model, or standard ACDM model Q. 
It is characterized by an essentially zero value of the spatial curvature parameter and a 
non- vanishing, but positive, value of the cosmological term (CT), A, in Einstein's equa- 
tions. The ACDM model supersedes the old Einstein-de Sitter (or critical-density) model, 
which has remained as the standard cold dark matter (CDM) cosmological model for many 
years. The CDM model has also zero curvature but at the same time it has zero CT 
(A = 0), in contradiction with present observations. Evidence for the values of the cos- 
mological parameters within the new concordance model comes both from tracing the rate 
of expansion of the Universe with high-z Type la supernovae experiments and from the 
precise measurement of the anisotropies in the cosmic microwave background (CMB) radi- 
ation pL t| ffl. These experiments seem to agree that 70% (0° ) of the total energy density 
of the Universe at the present time should be encoded in the form of vacuum energy 
density p\ or, in general, as a new form of energy density, po, which has been dubbed 
"dark energy" (DE) because of our ignorance about its ultimate nature H . The remaining 
30% {Q° M = ^cdm + ^b) °f the- present energy budget is essentially made out of cold dark 
matter, with only a scanty Q^ < 5% consisting of baryons. 

At present the ACDM model is perfectly consistent with all existing data [||, y, |ij , 
including the lately released three year WMAP results ||. In particular, this means that a 
non- vanishing and constant value A = Ao does fit the data reasonably well |Hj. However, 
when analyzing the pressure and density (pu , pd) of the mysterious DE fluid (presumably 
responsible for the accelerated expansion of the Universe), some recent studies jn], [l2|, [l3| 
seem to suggest that its equation-of-state (EOS) parameter, up = Pd/pd, could be evolving 
with time or (cosmological) redshift, uj£> = lud(z), showing a significant departure from the 
naive expectation lod = — 1 associated to the CT. Although this potential variation seems 
to be dependent on the kind of data set used, at present it cannot be ruled out [|14|, H. In 
some cases it could even accommodate lod < — 1, namely a "phantom- like behavior" (see 
below) jnj. In the near future we expect to get more accurate measurements of the dark 
energy EOS from various independent sources, including galaxy clusters, weak gravitational 
lensing, supernova distances and CMB anisotropies at small angular resolution. The main 
experiments that should provide this information are SNAP |[ , PLANCK || and DES || . 
If the overall result from these experiments would be that the DE density is constant 
(i.e. not evolving at all with the cosmological redshift), then p£> could just be identified 
with the constant energy density p\ ~ 10 -47 GeV 4 associated to the cosmological constant 
Ao = 8 ir G p\ of the standard ACDM model. Even though this could be looked upon as the 
simplest and most economical explanation for the DE, it should not be regarded as more 
satisfactory, for it would not shed any light on the physical interpretation of A in Einstein's 
equations. In particular, the notion of A as being a constant vacuum energy density 
throughout the entire history of the Universe clashes violently with all known predictions 
from quantum field theory (QFT), including our cherished Standard Model of strong and 
electroweak interactions. This leads to the famous cosmological constant problem (CCP), 



the biggest conundrum in Theoretical Physics ever [16, 17, 18, 19 1 . While we do not have 
at present a clear idea of how to associate the value of the vacuum energy to the CT, a 
dynamical picture for the DE looks much more promising. For, if pp = pd(S>) is a function 
of the cosmic time, or equivalently of the cosmological redshift pa = pu(z), then there 
is a hope that we can explain why it has the particular value p® D ~ 10~ 47 GeV 4 we have 
measured today and why its value may have been very different from the one (possibly 
much larger) that pjy took at early times, and the (tinier?) value it may take in the 
remote future (including a potential change of sign). This possibility has been suggested in 
different contexts [E3, E3] as a means for alleviating certain problems with the formulation 



of asymptotic states in string theory with positive A [24]. 

The notion of a variable vacuum energy is rather old and it was originally implemented 



in terms of dynamical scalar fields [25]. For example, the "cosmon" field introduced in [26] 
aimed at a Peccei-Quinn-like adjustment mechanism based on a dynamical selection of the 
vacuum state at zero VEV of the potential, < V >= 0. More recently these ideas have 
been exploited profusely in various forms, such as the so-called "quintessence" scalar fields 



and the like [27, 28], "phantom" fields [15| , braneworld models |2JJ, Chaplygin gas [30], and 
many other recent ideas like holographic dark energy, cosmic strings, domain walls etc 



(see e.g. [17, [Tg, [L9| and references therein), including some recently resurrected old ideas 



on adjusting mechanisms |pl] 1, and also possible connections between the DE and neutrino 



physics [32, 34] or the existence of extremely light quanta [32, 33 1. If the equation of state 
Pd = ujd pd describes some scalar field x (perhaps the most paradigmatic scenario for the 
dynamical DE models), and the EOS parameter up lies in the interval —1 < ujd < —1/3, 



then the field \ corresponds to standard quintessence [28|; if, however, lod < —1 then x is 
called a "phantom field" [[0| because this possibility is non-canonical in QFT (namely it 
enforces the coefficient in its kinetic energy term to be negative) and violates the weak and 



(o fortiori) the dominant energy conditions [35]. 

Let us emphasize that the EOS parameter u>d could be an effective one |3g, g7|> g8[ , 
and hence we had better call it u e . Actually, the effective uj e could result from a model 
which does not even contain a single scalar field as the direct physical support of the DE. 
For example, in Ref.|37]] it has been shown that a running A model leads to an effective 
EOS, po = u e p£j, which may behave as quintessence or even as a phantom-like fluid. 



Remarkably this feature has been elevated to the category of a general theorem [38], which 
states the following: any model based on Einstein's equations with a variable A = Kit) 
and/or G = G(t) leads to an effective EOS parameter u) e which can emulate the dynamical 
behavior of a scalar field both in "quintessence phase" (— 1 < uj e < —1/3) or in "phantom 
phase" (u e < — 1) - therefore always entailing a crossing of the cosmological constant divide 
— 1 |pq| . For another example of this general theorem, in a complementary situation 



U). 



where G is variable but p\ is constant, see |39|| . 

In spite of the virtues of cosmologies based on variable cosmological parameters, it is 
convenient to study the possibility of having a composite DE model involving additional in- 
gredients. As we shall see, this may help to smooth out other acute cosmological problems. 



1 See also Ref. J20L Elf for a summarized presentation of the CCP. 



For instance, consider the so-called "cosmological coincidence problem" [4C], to wit: why do 
we find ourselves in an epoch t = to where the DE density is similar to the matter density 
(PD(to) — /9m(*o))? I n the ACDM model this is an especially troublesome problem because 
p\ remains constant throughout the entire history of the Universe. In the ordinary dynam- 
ical DE models the problem has also its own difficulties. Thus, in a typical quintessence 
model the matter-radiation energy density p m = pu + PR decreases (with the expansion) 
faster than the DE density and we expect that in the early epochs p m ^> po, whereas at 
present and in the future p m <C pd- Therefore, why do we just happen to live in an epoch 
to where the two functions /3d (to) — Pm(^o)? Is this a mere coincidence or there is some 
other, more convincing, reason? The next question of course is: can we devise a model 
where the ratio p m /pD stays bounded (perhaps even not too far from 1) in essentially the 
entire span of the Universe lifetime? This is certainly not possible neither in the standard 
ACDM model, nor in standard quintessence models 2 . And it is also impossible for a model 
whose DE consists only of a running A [p3fl . However, in this paper we will show that we 
can produce a dynamical DE model with such a property. Specifically, we investigate a 
minimal realization of a composite DE model made out of just two components: a running 
A and another entity, X, which interacts with A. In this framework matter and radiation 
will be canonically conserved, and the running of A is (as any parameter in QFT) tied to 
the renormalization group (RG) in curved space-time [g(], |44J| . To illustrate this possibility 
we adapt a type of cosmological RG model with running A that has been thoroughly stud- 
ied in the literature [£J| 3 . On the other hand a most common possibility for X in QFT and 
string theory would be some scalar field \ ( e -g- moduli or dilaton fields |£^]) that results 
from low-energy string theory. In fact, the original cosmon field was a pseudo-dilaton p6|j . 
Here we will also use the "cosmon" denomination for X. This entity will stand for dynam- 
ical contributions that go into the DE other than the vacuum energy effects encoded in A. 
Since some of the effects associated to X could be linked to dynamical scalar fields, this 
justifies its denomination of cosmon [26|. Whatever it be its ultimate nature, we require 



that the total DE density pd must be conserved. We will show that the resulting AXCDM 
model can provide a good qualitative description of the present data, particularly the fitted 
EOS, and at the same time it helps to tackle the cosmological coincidence problem. 

The structure of the paper is as follows. In the next section we describe a few general 
features of composite DE models and construct the AXCDM model. In Section || we explic- 
itly solve the model. In Section ^| we consider the constraints imposed by nucleosynthesis 
and the evolution of the dark energy. In Section |E| we describe possible cosmological scena- 
rios within the AXCDM model. In Sect. |6| we explore the possibility of crossing the uj e = — 1 
divide. In Sect. fj\ we perform a comprehensive numerical analysis of the AXCDM model 
and discuss how this model can help to alleviate the cosmological coincidence problem. In 
the last section we present the final discussion and we draw our conclusions. 



2 There is, however, the possibility to use more complicated models with interactive f4T| or oscillating |42] 
quintessence (see also [|19|), or to consider the probability that pn/ 'pm ~ 1 in a phantom Universe [ |43| . 

^For the recent literature on RG models supporting the idea of running cosmological parameters, whether 
in QFT in curved space-time or in quantum gravity, see e.g. [4q, ha, WA km . See also M, pOI] for various 
phenomenological models with variable cosmological parameters. 



2. AXCDM model: dynamical DE made out of A and cosmon. 

The simplest notion of dynamical dark energy can be implemented in terms of a single 
dynamical field X that completely supersedes the cosmological constant. This point of 
view gave rise to the notion of XCDM model M as a substitute for the standard ACDM 
one [y . When considering a XCDM model with a single DE fluid of general kind (even 
a phantom one), and under the general assumption < O^ = p° M / ' p° c < 1, the energy 
density po must be positive at present, p D > 0, otherwise it would be impossible to fulfill 
the cosmic sum rule for flat universes: Q? M + Q D = 1, with Q D = p D /Pc- Furthermore, 
one usually assumes that this situation is always the case, i.e. po > at any time in 
the cosmic evolution. This entails that only the DE pressure, pa, can be negative. For 
example, for a quintessence field we have pn > and —po < Pd < — Pd/3, so that the 
strong energy condition (SEC) is violated, but not the weak (WEC) and dominant (DEC) 
ones (cf. Fig.[j]). However, it should be noted that the data is usually described by an 
overall effective EOS, po = ^ePD, irrespective of the potential composite nature of the 
fluids that build up the total DE density. In composite DE models, with a mixture of 
fluids of individual EOS's pi = uJi pi (i = 1, 2, ..., n), the effective EOS of the mixture reads 

Pd wipi +u 2 P2 + ■■■ , 01 v 

u e = — = . (2.1) 

PD P1+P2 + --- 

Clearly, uj e is not related in a simple way to the various Wj and in general to e = co e {z) will be a 
complicated function of time or redshift even if all tOi are constant. A mixture of barotropic 
fluids is, therefore, non-barotropic in general. As a consequence many possibilities open 
up. For example, it should be perfectly possible in principle to have a mixture of fluids 
in which one or more of them have negative energy density, pi < 0, and positive pressure 
Pi > 0, even at the present time. These DE components can be phantom-like, in the sense 
that may result in Wj < —1, but they are not phantoms of the "standard" type, for they 
violate the WEC and DEC while they can still preserve the SEC. Such "singular" phantom 
components are very peculiar because they actually reinforce the cosmological gravitational 
pull, in addition to the one exerted by ordinary matter (see Eq. ( |2.13| ) below). In this 
sense, rather than acting as conventional phantom energy, these components behave as a 
sort of unclustered "matter" (with negative energy) that may be called "Phantom Matter" 
(PM) (cf. Fig.|l|). Although the latter is a bit bizarre, the rationale for it will appear 
later. The point is that these components may eventually cause the halt of the expansion. 
Phantom Matter could be an unexpected component of DE, and one that still preserves 
the SEC as normal matter does. In contrast, a standard phantom energy component 
realizes the condition u>i < —1 in the form pi < —pi < (see Fig. Q). It violates all 
classical energy conditions (WEC, DEC, SEC, including the null and null-dominant energy 
conditions p3]) and produces a specially acute anti-gravitational effect (much stronger 
than that of a positive cosmological term) which leads eventually to the so-called "Big 
Rip" singularity, namely the Universe is ripped to pieces after a super-accelerated phase 
of finite duration |ll| . Therefore, the two forms of phantom-like behavior o; e < — 1 are 
dramatically different: whereas standard phantom energy produces super-acceleration of 
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Figure 1: Energy conditions: (a) The shaded regions fulfill the Weak Energy Condition (WEC): 
p > and p + p > 0. The lighter-shaded one also satisfies the Dominant Energy Condition (DEC): 
P > \p\] (b) The region shaded in gray (with or without dots) fulfills the Strong Energy Condition 
(SEC): p + p > and p + 3p > 0. The quintessence (Q) region (—1 < tu e < —1/3) is marked cross- 
hatched, and the usual phantom region (uj e < — 1 with p > 0) is indicated by P. The gray-dotted 
region corresponds to an unusual kind of phantom uj e < — 1 with p < 0, which we call "Phantom 
Matter" (PM). Note that PM satisfies the SEC. The cosmon behaves in many cases as PM within 
the AXCDM model. 



the expansion and leads to Big Rip destruction of the Universe, Phantom Matter produces 
super-deceleration and fast stopping (with subsequent reversal) of the expansion. That 
said, it is obvious that what really matters in composite DE models is the behavior of the 
overall parameters (pd,Pd)- Irrespective of the particular nature of the DE components 
we can still have e.g. the overall density pn positive and the total effective pp negative. 
Then if — 1 < tu e < —1/3 the effective behavior of the DE fluid mixture is quintessence-like, 
while if <jj e < — 1 it is phantom-like (of the standard type), such that the cosmic sum rule 
for flat universes can be preserved. 

As we have mentioned, experimentally an effective (standard) phantom fluid (po > 0, 
uj e < — 1) for the DE cannot be discarded at present because it is allowed by the combined 
analysis of the supernovae and CMB data 0, O, [13]. In the following we will construct a 
simple composite model where this kind of mixture situations can be met. The resulting 
AXCDM model is a kind of "composite XCDM model" because, in contradistinction to the 
original XCDM ||, we keep both the cosmological term (with a given status of dynamical 
entity) and also the new dynamical fluid X, and we assume that in general they are in 
interaction 4 . Obviously a negative value of the cosmological term at present (Ao < 0) is 
allowed in this framework (with the current total DE density po > 0) thanks to the X 
component. We will see that in the future this may cause a halt of the expansion. This 
issue has been investigated also in models with a fixed negative value of the cosmological 
term p2]. However, thanks to the variability of A and/or the PM behavior of X, the halt 



4 One might entertain whether one should better call it XACDM, rather than AXCDM. In the former 
case we would seem to imply that we introduce X to the ACDM model |l[ , whereas in the latter we would 
be "restoring" A into the XCDM one M. However, we emphasize that the presented model is not just the 
result of addition, but also of interaction, and in this sense the name should not be a problem. We stick to 
AXCDM, with no subtle connotations whatsoever. 



will also be possible for Ao > 0. The combined dynamics of A and cosmon X opens many 
other possibilities that cannot be offered by models which just add up a quintessence field 
to the standard ACDM and commit to some particular form of the potential for that field. 
The fact that we do not specify the nature of X by e.g. identifying it with a scalar field, 
endows our discussion with a higher degree of generality. 

It should be stressed right from the start that despite the time variation of the cos- 
mological parameters within our framework, this fact in no way contradicts the general 
covariance of the theory. This is independent of whether such variation is inherited from 



RG considerations or it has some general unspecified origin [49, 5Cj. The parameters could 
even be variable in space as well (see e.g. [flq]), although we will not allow them to do so 
here because we want to preserve the Cosmological Principle. The general covariance is in- 
sured by imposing the fulfillment of the general Bianchi identity on both sides of Einstein's 
equations. The kind of particular relations that this identity will imply among the various 
cosmological quantities will depend on the metric itself, which in our case we choose to 
be the Friedmann-Lemaitre- Robertson- Walker (FLRW) one 5 . The most general form of 
the Bianchi identity expressing the equation of continuity for a mixture of fluids (^T|) and 
variable cosmological parameters reads: 



d 
~dl 



g(Ep> 



+ GH^2a iPl = 0, a t = 3(l+iOi). (2.2) 



In our case we assume that we have matter-radiation, p m = pu + PR, and two DE compo- 
nents: one is called X ("cosmon"), with dynamical density px = Px(t), and the other is a 
cosmological term with density p\, which can be constant but in general it is also allowed 
to vary with the evolution: p\ = ph{t)- This is perfectly possible |Q] and preserves the Cos- 
mological Principle provided it only varies with the cosmic time or (equivalently) with the 
cosmological redshift: p\ = p\{z). The two DE components have barotropic indices ujx 
and uj\ = — 1 respectively. In general the index ux for the cosmon will be dynamical, but 
in practice one can assume that X is some barotropic fluid with constant lux, typically in 
one of the two expected ranges wj > — 1 (quintessence- like) and wj < — 1 (phantom- like) . 
For the moment we assume it is an arbitrary function of time or redshift. It should be 
clear that for the "A fluid" the relation uj\ = — 1 holds irrespective of whether /?a is strictly 
constant or variable. In this work we will keep the gravitational coupling G constant and 
assume that the matter-radiation density p m = pu + PR is conserved. Hence Eq. ( |2.2j ) 
splits into corresponding conservation laws for matter-radiation and total DE, namely 

Pm+a m p m H = 0, a m = 3(1 + u m ) , (2.3) 



J A different and certainly independent issue is whether these models admit in general a Lagrangian 
formulation, and in this sense they are phenomenological insofar as the connection to that formulation is 
not straightforward. But in our opinion this should not necessarily be considered a priori as a shortcoming. 
After all the original introduction of the field equations (by Einstein himself) was based purely on general 
covariance considerations and with no reference whatsoever to an action principle. This said, a Lagrangian 
formulation of some of these models is not excluded and may lead to different phenomenological conse- 
quences. From our point of view this issue must be decided by experiment. For instance, in Ref. p3| it was 
shown in great detail how to test one of these RG models using the distant supernovae data. 



with 

Um = \ P l , (2.4) 

3 PM + PR 

and 

p D + a e p D H = 0, a e = 3(l+uj e ), (2.5) 

with 

PD = PA + px- (2.6) 



In Eq. (|2.4| ) w m = 0,1/3 (a m = 3,4) for the matter dominated epoch and the radiation 
dominated epoch respectively. Clearly, the existence of a dynamical DE component X is a 
necessary condition for p\ being variable in this framework. In the particular case where 
PA is constant, X would be self-conserved, but in general pa and X may interact and 
therefore they can constitute two different dynamical components of the total DE density 
( |2.6| ). The effective barotropic index - or effective EOS parameter fl2.1|) - of the AXCDM 

model is 

PA+Px -PA + ^xPx , . ,-. . ,px ,„ ^ 

ui e = ■ = ■ = -1 + (1 + ujx) — • (2.7) 

PA + PX PA + PX PD 

As noticed above, even if X is a barotropic fluid with constant wj the mixture of X and 
A is non-barotropic because px and /?a will in general be functions of time or redshift, and 
so will be uj e = oj e (z). We shall compute this function precisely for the AXCDM model 
under certain dynamical assumptions on the behavior of the cosmological term. From the 
equations above it is easy to see that the overall DE conservation law ( |2.5D can be rewritten 
as follows: 

PA + Px + a x pxH = 0, a x = 3(l + uJx)- (2.8) 

For constant pa, this law boils down to the self-conservation of the X component, 

px+a x pxH = 0. (2.9) 

For variable pa, Eq. (|2.8|) shows that the dynamics of the cosmological term and the cosmon 



become entangled. One can confirm Eq. (2.8) starting from the total energy-momentum 
tensor of the mixed DE fluid, with 4- velocity U^: 

T° = rf v + T* v = (p A - wi p x ) 9»» + (1 + lo x )px U„ U v . (2.10) 

Then one can use the FLRW metric 

ds 2 = dt 2 _ 2( t ) f -^ + r 2 d6 2 + r 2 sin 2 9 dA (2.11) 

\1 — kr z J 

and straightforwardly compute \J^T D ^ = 0. The result of course is fl2.8|) . Another fun- 
damental equation of the cosmological model under construction is Friedmann's equation 

„•> { a\ 8ttG, . k 8ttG, . k ,„ .,„, 

H =(-) =^-(p m + PD)- — = ^-(pm + PA + Px)- — - (2.12) 

\ (X / o do (X 



Let us also quote the dynamical field equation for the scale factor: 

a 4:ir G , . . . . 47r G , 

- = 5— [Pm (1 + 3iO m ) + p D (1 + 3w e )J = — [p m (1 + 3iO m ) - 2 p A + px (1 + 3^x)J • 

(2.13) 

This one is not independent from (|2.12 ) and ( |2.8[) , as can be easily seen by computing the 

time derivative of H with the help of equations ( |2.12| ) , (£l]) and ( |2.3| ) : 

_ v 47rG, , k 4-kG . , fe .„ n ^ 

-n = —[a m Pm + a e p D \ + ^j = — a m p m + axPx + -0 • ( 2 - 14 ) 

3 a z 3 or 

Then Eq.( 2.13| ) immediately follows from a/ a = H + i7 2 , showing that indeed it is not 
independent. However it is useful to have ( 2,13| ) at hand because it exhibits the acceleration 
law for the scale factor under the combined influence of the matter-radiation density p m 
and the composite DE density (2.6). 

How do we proceed next? Clearly, even if the index ujx for the cosmon would be 
given, we have three independent equations, e.g. (|2.3D , (|2.8|) and Q2.12 ), but four unknown 
functions (H (t) , p m (t) , p\(t) , px (t)) . As a matter of fact, Eq.(^) is decoupled from the 
rest, and its solution in the matter dominated and radiation dominated eras can be cast 
in the simple unified form 

p m (z) = p° m (l + z) a ™ , (2.15) 

where a m = 3,4 respectively for each era. For convenience we have expressed this solution 
in terms of the cosmological redshift, z, with the help of H dt = —dzj(\ + z). Notice that 
Eg. ([2.5]) can also be solved explicitly for the total DE density: 

p D (z) = pi exp J3 jT rf ^+^f) } . (2.16) 

This equation is basically formal. In practice it cannot be used unless the model is explicitly 
solved, meaning that we need first to find the individual functions p\ = p\{z) and px = 
px(z) for the two components of the DE. In particular, this equation does not imply that 
Pd = Pd(z) is positive definite for all z because, after solving the model, the function uj e = 
u) e {z) may present singularities at certain values of z. The reciprocal relation expressing 
the effective EOS parameter in terms of the total DE density is 

l \ 1 i 1 + Z l dpD in 1 v\ 

u e (z) = -l + — t-. (2.17) 

3 pd dz 



From (U) and (|2.17|) we get 

a x px= ( 1 + Z )^T ( 2 - 18 ) 

which shows that the effective quintessence {dp^/dz > 0) or phantom-like (dp^jdz < 0) 
character of the AXCDM model will be known once the sign of ax and the sign of the 
density px of the cosmon field are known. In particular, we note the curious fact that even 
if ujx < — 1 the model could still effectively look as quintessence-like provided we admit the 
possibility that the cosmon field can have negative energy-density {px < 0) 6 . In general 

6 The possibility of having negative energy components has been entertained previously in the literature 
within e.g. the context of constructing asymptotically de Sitter phantom cosmologies |52]. 



the AXCDM model effectively behaves as quintessence (resp. phantom DE) if and only if 
ax and px have the same (resp. opposite) signs. 

In the end we see that we have two coupled equations ( |2.8| ) and ( |2.12 ) for three un- 



knowns (H(t), p\(t), px(t)). We need a third equation to solve the AXCDM model. Here 
we could provide either a specific model for X or one for p\. Since X represents a generic 
form of dynamics other than the cosmological term (the scalar fields being just one among 
many possible options), from our point of view it is more fundamental to provide a model 
for pa. The simplest possibility is to assume p\ = const. In this case the new model is a 
trivial extension of the standard ACDM one, but even in this case it leads to interesting 
new physics (see scenario I, Section ||). However, in the general case pa could be a dynami- 
cal quantity and a most appealing possibility is to link its variability to the renormalization 
group (RG) [|(]] . Once the general principles of QFT in curved space-time (or in Quantum 
Gravity [^]) will determine the RG scaling law for A, the dynamics of X will be tied to the 
conservation law (^^) in a completely model-independent way; that is to say, independent 
of the ultimate nature of the cosmon entity. This will be true provided we adhere to the 
commonly accepted Ansatz that the total DE must be a conserved quantity independent of 
matter. While we do not know for the moment what are the general principles of QFT in 
connection to the running properties of the cosmological term in Einstein's equations, we 
have already some hints which may be put into test. As a guide we will use the cosmologi- 



cal RG model described in [23], based on the framework of [20, 46 1 within QFT in curved 
space-time. A more general RG cosmological model with both running G and running 
PA can also be constructed within QFT in curved space-time, see Ref . |45f| . However, for 
simplicity hereafter we limit ourselves to the case G = const. Briefly, the model we will use 
is based on a RG equation for the CT of the general form 

M^^^' (2 ' 19) 

n=l 

Here p is the energy scale associated to the RG running in cosmology. It has been argued 
that p can be identified with the Hubble parameter p = H at any given epoch p0|. Since 
H evolves with the cosmic time, the cosmological term inherits a time-dependence through 
its primary scale evolution with the renormalization scale p. Coefficients A n are obtained 
after summing over the loop contributions of fields of different masses Mj and spins U{. 
The general behavior is A n ~ ^2 M t ~ n [20, 46]. Therefore, for p <C Mj, the series above 
is an expansion in powers of the small quantities p/Mi. Given that A\ ~ YL^ii the 
heaviest fields give the dominant contribution. This trait ("soft-decoupling") represents a 
generalization of the decoupling theorem in QFT- see [p3J | for a more detailed discussion. 
Now, since p = Hq ~ 10 -33 eV the condition p <C Mj is amply met for all known particles, 
and the series on the r.h.s of Eq. (2.1S|) converges extremely fast. Only even powers of 



p = H are consistent with general covariance [23|, barring potential bulk viscosity effects 
in the cosmic fluid [|3J|] . The n = contribution is absent because it corresponds to terms 
oc Mf that give an extremely fast evolution. Actually from the RG point of view they are 
excluded because, as noted above, p <C Mj for all known masses. In practice only the first 
term n = 1 is needed, with Mj of the order of the highest mass available. We may assume 



that the dominant masses Mi are of order of a high mass scale near the Planck mass Mp. 
Let us define an effective mass M as the total mass of the heavy particles contributing to 



the /3-function on the r.h.s. of (2.19) after taking into account their multiplicities aj [23|: 



M = 



\ 



E a * M ? 



(2.20) 



We introduce (as in J23J) the ratio 



a M 2 



(2.21) 



Here a = ±1 depending on whether bosons or fermions dominate in their loop contributions 

to ( p. 19 ). Thus, if the effective mass M of the heavy particles is just Mp, then v takes the 

canonical value avo, with 

1 

2.6 xl(T 2 . (2.22) 



i/ = 



12 7T 



We cannot exclude a priori that v can take values above ( 2.22J ) if these multiplicities are 
very high. Under the very good approximation n = 1, and with \i = H, Eq. ( [2.19D abridges 
to ||| 

^ = ^M 2 P H 2 . (2.23) 

d\nH An ' 

Equation ( 2.23| ) is the sought-for third equation needed to solve our cosmological model. 

Its solution is 



with 



co 



PA 



PA 

3v 

8^ 



c + ci H , 



M P Hi , d 



3^ 



■M 



-in 



p ■ 



(2.24) 



(2.25) 



where H = H(t) is given by Q2,12j ). Notice that this result gives some theoretical basis to 
the original investigation of the possibility that A ~ H 2 (equivalently, pa ~ H 2 Mp) |l|] at 
the present epoch. However, in our case rather than a law of the type p\ ~ H 2 Mp we 
have 5 p\ ~ H 2 Mp, and the variation 5p\ around the reference value is controlled by a 
RG equation. 



3. Analytical solution of the AXCDM model 



Using ( |2.8D ( 2.12p , ( p. 14 ) and ( 2.23J ) we can present the relevant set of equations of the 
AXCDM model in terms of the cosmological redshift variable z as follows: 



dp x dp A 



ax Px 



dz 

dpA 

dz 

dH 2 

dz 



dz 

3u 



1 + 2 

. ,2 dH 2 
Mp 



8vr ^ dz ' 
- 8 -^ amPm + axPX +2H 2 n° K (l + z) 



(3.1) 
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The latter can be replaced by Friedmann's equation ( |2.12| ) at convenience. Here f^ = 
—k/H^ is the space curvature term at the present time. 

Equations ( |3.1|) must be solved for (H, px, Pa), where p m is known from Eq. ( 2.15| ). 
Instead of p\ , px we will use 



Q A (z) 



Pa(z) 



n x 



Px(z) 

P°c 



where p a c = 3Hq/87tG is the critical density at present. We also define 



^■' j m\Zj 



Pm(z) 



PY 



n K (z) = n° K (i + z) 



From these equations we can derive the following differential equation for Six'- 



dsi 



x 



dz 



(l-i/) 



a x Six 



Oim, »"rr 



2un° K (i + z ) 



(3.2) 



(3.3) 



(3.4) 



1+2 l+Z 

We can solve it for Six exactly, even if ujx (hence ax) is an arbitrary function of the 
redshift. The result is the following: 



Slx(z) ~- 


-a\ 


1- / dz" q 

Jo 


{z") exp 




l + Z 






q(z) = v 


r n° 



dz p(z ) 



cxp 



dz p{z ] 



(3.5) 



1 + z 



iOm-1 



n 



2 4 (1+ - 



(3.6) 



Notice that if v = then Sl\ is strictly constant, Eq. (|3.4j) reduces to Eq. ( ^.9[ ) and Slx(z) 
becomes self-conserved, as expected. If, in addition, ujx =const. then (|3.5|) collapses at 
once to 



Slx{z) = Sl x {l + z, 



ax 



(3.7) 



which is of course the solution of ( |2.9| ) for constant ax- This situation is the simplest pos- 
sibility, and we shall further discuss it in Section |5| as a particular case of our more general 
framework in which the two components of the DE are dynamical and interacting with 
one another. The next to simplest situation appears when v ^ (still with lux =const.) 
because then there is an interaction between the two components of the DE and these 
components can be treated as barotropic fluids. In such case we can work out the integrals 



in (|3.5| ) and obtain a non-trivial fully analytical solution of the system (34) for arbitrary 
spatial curvature. The final result can be cast as follows: 

HHz) ... . 



Hi 



Sl A (z) 





- V 


1- 


V 


Ji A 


- V 



l-v 

Slx(z) = {l-v)F{z) 



+ uF{z), 



ST m (l + z) 



Sl° K (l + z) 2 . 



(3.8) 
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with 



F(z) 



n° m (a r 



ax) 



:i+z) 



a, 



+ 



ax (1 - v) 



"8r(°x-2) , 



ax (1 - ^) 



i-o°a n°ic. 



ax) , ^(«i - 2) 



1-1/ 



+ 



a, 



ax (1 - i/) 2 - ax (1 - ^) 



(l + ? )ojf(H. (3.9) 



The solution above corresponds to the matter (a m = 3) or radiation (a m = 4) dominated 
eras. The expansion rate for the AXCDM model in any of these eras can be expressed as 



H\z) = Hl [Q^(l + z) am +n° K (l + z) 2 + n D (z)] , 



(3.10) 



where Qjj(z) = pd(z)/p® is the total DE energy density in units of the current critical 
density. From the previous formulae the latter can be written in full as follows: 



n D (z) 



o° 



V V 



a x n° m (l + z) a ™ ua x n° K (l + zf 



+ 



/' a r 

)° 



a x (1 - v) 

L m\ a m 



2-ax{l-y) 
l-to\ _ f4(q m - a x ) ^ _V° K (ax - 2) 
1 — v a r 



+ 



ocx (1 - v) 2 - a x (1 - v) 



(1 + *) 



a x (1— v) 



. (3.11) 



Let us finally derive the generalized form of the "cosmic sum rule" within the AXCDM 
model. From ( p.ll ) we can readily verify that for z = we have ^^(O) = 1 — ^m~ ^k- This 
renders the desired sum rule satisfied by the current values of the cosmological parameters 
in the AXCDM: 



M ' D ' 1 






Q° A + 0° + 



The cosmic sum rule can be extended for any redshift 

n m ( z ) + n A (z) + n x (z) + n K (z) = 1 

provided we define the new cosmological functions 

Px(z) 



O / X P™{z) A . v gA(g) o / . 

il m (z) = — — , Ua(z) = — — , ilx{z) 



Pc{z) 



Pc(z) 



Pc{z) 



, £Ik{z) 



1. (3.12) 



(3.13) 



(1 + z) 2 , (3.14) 



# 2 (z) 



with p c (z) = 3H 2 (z)/8ir G the critical density at redshift z. These functions should not be 
confused with ( |3.2[ ) and ( |3.3| ). We will, however, use more frequently the latter. 

We may ask ourselves: in which limit do we recover the standard ACDM model? We 
first notice that for v = the function fl3.9|) dwarfs to 



f(z) = n° m (i + zY™ + n° K (i + z) 2 + o x (i + *; 



"X 



(3.15) 



where we have used ( 3.12 ). Then substituting this in ( |3.g| ) one can immediately check that 
the standard ACDM model formulae are obtained from those of the AXCDM model in the 
limit v — > and fi^ — > 0, as could be expected. 
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4. Nucleosynthesis constraints and the evolution of the DE from the ra- 
diation era into the matter dominated epoch 

At temperatures near 1 MeV the weak interactions (responsible for neutrons and protons 
to be in equilibrium) freeze-out. The expansion rate is sensitive to the amount of DE, 
and therefore primordial nucleosynthesis (T ~ 0.1 — lMeV) can place stringent bounds 
on the parameters of the AXCDM model. A similar situation was considered in the run- 
ning cosmological models previously studied in Ref. [^, [37J| . In the latter, nucleosynthesis 
furnished a bound on the parameter v defined in Eq.( |2,2l| ) from the fact that the ratio of 
the cosmological constant density to the matter density should be relatively small at the 
nucleosynthesis epoch, (p\/p m )N <C 1, of order 10% (see also (54]]). This condition led to 
v < 0.1, which is e.g. amply satisfied by the natural value (|2.22|). However, in the AXCDM 



model the total DE density is not just p\, but pr> given by (2.6), and moreover the dynam- 
ics of the cosmological term here is not linked to the matter density (as in Ref.[23|) but to 
the cosmon density. Therefore, the appropriate ratio to be defined in the present instance 

PD PA + PX "A + "X ,. 1N 

r = — = — " = — o • I 4 - 1 ) 

Pm Prn * '"m 

Notice that for the flat case {k = 0) this ratio is related to Clo = ^A + &x (cf. Eq.( 3.14| )) 
as follows: 

n D = JL- . (4.2) 

Even if k ^ 0, at the nucleosynthesis time the curvature term is negligible. Therefore 
the nucleosynthesis bounds on Qd of order 10% 54] are essentially bounds on r. Let us 



compute the ratio r in general and then evaluate it at nucleosynthesis. For convenience let 
us use the function (|3.9[) and express the result as 



r ( z ) = on r, , .a»„ • ( 4 - 3 ) 



+ F{z) - UP m (1 + z) a ™ -n u K {l + z 

nui + z) am 

It is easy to check that for z = this equation reduces to the expected result 



r = r ( 0) = ^i = ^. (4.4) 

Let now z N ~ T N /T ~ (0.1 MeV)/ (2 x 10~ 4 eV) ~ 10 9 be the typical cosmological redshift 
at nucleosynthesis time. Let also ojr = 1/3 and aji = 4 be the values of the parameters 
uj m and a m at nucleosynthesis (and in general in the radiation epoch). We can obviously 
neglect the first (constant) term in the numerator of ( [4.3[ ) at z = z^r, and as we have said 
also the contribution from the spatial curvature (which scales with two powers less of (1 + z) 
as compared to the radiation density). After a straightforward calculation we obtain 

F(z N )-n° R (i + ZN ) A 

r N = r(z N ) ~ f =r e + R N (z N ) (4.5) 

where 

r e = '——, (4.6) 

uj R -uJx + e 
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with 



and we have defined the function 



e = v(l + wx) 



(4.7) 



Rn{zn) 



n% 



wb -«x 



Q%(1 - v) lo r -uj x + e 



<& 

^ 



1 + 3wx 
1 + 3 wx — 3 e 



(I + ztv) 3 ^"^" 1 - (4. 



Here Q R is the radiation density at present. The parameter e defined in (4.7) will soon play 
an important role in our discussion. The "residual" function ( f4.8[ ) will be totally harmless 
at nucleosynthesis provided the exponent of (1 + zx) that goes with it is sufficiently large 
(in absolute value) and negative, in which case Rn{zn ~ 10 9 ) will be negligible. We start 
with the condition that it be negative, which provides the constraint 



LU X 



1 

£< 3- 



(4.9) 



However, there is also the first term on the r.h.s. of ( |4.5| ). This term does not decay 
with the redshift, it is constant and therefore must be small by itself. If we take the same 
condition as in [23] concerning the bound on the ratio of DE density to radiation density 
at nucleosynthesis, namely |rjy| < 10%, we find the additional constraint 



(4.10) 



<0.1 . 

ur - lux + e 

We shall assume that the cosmon barotropic index lies in the interval 

-1 - 5/3 <uj x < -1/3 (5 > 0) , 



(4.11) 



with 5 a small positive quantity. The upper limit is fixed by requiring that X has at 
most quintessence behavior (i.e. excluding matter-radiation behavior); the lower limit is 
expected to be some number below (but around) —1 in order to accommodate a possible 
phantom-like behavior. From Eq. ( |4.11[ ), and recalling that ujr = 1/3, it is patent that the 
bound ( 4.1CJ ) just translates into a rough bound on e itself: 



\r e \ — e 



\v(l+LU X )\ < 0.1. 



(4.12) 



Using this relation and the condition (4.11) we see that the previous bound fl4.9|) is amply 



satisfied and the function (4.8) actually behaves as 

R N (z)~(l + z) (/3<-2-3e). 



(4.13) 



This means that Rn{zn ~ 10 9 ) will be completely negligible, as desired. Therefore, the 
two bounds from nucleosynthesis just reduce to one, Eq. ( 4.1C| ), and we can simply set 
there |rjv| = |?" e |. We will use the exact bound ( 4.10| ) for our numerical analysis, but the 
approximate version ( 4.12|) is useful for qualitative discussions. 



In contrast to the nucleosynthesis bound obtained in Ref. 23] for the alternative sce- 
nario where the DE was composed only of a running cosmological constant (in interaction 
with matter) we see that for the AXCDM model it is not necessary to require that \v\ is 
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small, but only the product of \v\ times \ax\- We also note that in the AXCDM model 
for v ^ there is an irreducible contribution of the DE to the total energy density in the 
radiation era, whose size is basically fixed by the tolerance of the nucleosynthesis processes 
to the presence of a certain amount of dark energy. The existence of this residual constant 
contribution shows that in the AXCDM model the presence of the DE component takes 
place at all epochs of the Universe evolution. It is therefore interesting to study what will 
be the evolution of this tiny (but not necessarily negligible) contribution in later times, in 
particular in the matter dominated epoch, where u) m = (a m = 3). Let us thus compute 
the ratio ( |P| ) in this epoch, and call it tm(z). In this case we apply the full formula (fO|), 
without neglecting any term, because z in the matter-dominated epoch covers a range from 
~ 10 3 to (near) — 1. The final result is: 



ru{z) 



n\ 



(i 



si Af 



with 



Rm(z) 



l + z) 



nl 



+ 






3e 



+ 



1 + 3o»x — 3e) (1 +z) ujx - e 



+ R M (z), (4.14) 



LU X 



n° M (i-u) ujx 



0° 

5i A7 



1 + 3 UJ X 

1 + 3 ujx — 3 e 



(1+^)3(^-6) 



It is easy to check that for z = the expression ( |4.14 ) shrinks just to 



(4.15) 



For general 



we see that the first term on the r.h.s. of ( 4.14| ) goes like (l + z) 3 and hence evolves 



fast with the expansion and is unbounded for z — ► — 1. The second term (the curvature 
dependent term) evolves much more slowly, namely as (1 + z)^ 1 and it can be neglected 
to study the future evolution (also because Cl^ i s ver y small). The rising of r(z) with the 
expansion cannot be stopped in the standard ACDM model because for it the ratio just 
reads 

r M(z) = on ,/, ,, ■ (4.16) 



n° M (i + 



This particular result is of course recovered from ( 4.14 ) in the limit v = and £l x = 
(see Section |||). However, in the AXCDM the presence of the non-trivial function Rm(z) ~ 
(1 + z) 3 ( w ^ _e ) is crucial as it allows the existence of both a stopping point, z s , as well as of 
an extremum of the function (4.14) at some z = z max during the future Universe evolution. 
We will study the relationship between these two features in the next section and we will 
see that the stopping condition is correlated with this extremum being a maximum. Here 
we will concentrate on the conditions of existence of this maximum of tm(z)- 

Let us from now on restrict our considerations to the flat space case, which is the most 
favored one. The redshift position z. 
The result is the following: 



of the maximum can be worked out from Eq. (4.14) 



l + z n 



O 



v 



K) 



3( 1+UJX - € ) 



(4.17) 



The nature of this extremum is obtained by working out the second derivative of r(z) at 
z = z max . We find 



V^ ~r ■2maxJ'~ 



0° 
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where z max is given by ( 4.1 7| ) . A sufficient condition for the extremum to be a maximum 
is that 

a x (tt° A -v) <0. (4.19) 

The height of the maximum, r max = r(z max ), is the following: 



n% 



V 1 +UJX 



(1 + z n 



0° 



+ 



LUX 



LU X 



(4.20) 



where the last term on the r.h.s. is negligible. 

Notice that the aforementioned maximum can only exist in the future (or at present), 
but not in the past. Its existence in the past would be incompatible with the experimentally 
confirmed state of acceleration of the Universe that started in relatively recent times. To 
prove this, one can show that the acceleration Eq. ( |2.13| ) can be written in terms of the ratio 
( |4.1| ) and its derivative r'(z) with respect to redshift as follows. We first write Friedmann's 
equation for k = in the form H 2 = (8ttG/3) (1 + r) p m . Then we differentiate it with 
respect to the cosmic time, 



H 



4vrG 



(1 + r) a m - — 



where use has been made of the equation of continuity (E 
equations we compute a /a = H + H 2 : 



(4.21) 



1). With the help of these 



a 

a 



AttG 

3 
4ttG 



(l + 3w m )(l+r(t)) 



r(t). 
H 



[(1 + 3w m )(l + r(z)) + (l + z) r'(z)] p n 



(4.22) 



where r'(z) = dr(z)/dz. The last equation shows that if the maximum could take place at 
some point z = z\ > in the recent past, then any z in the interval < z < z\ would satisfy 
r'(z) > 0, and then Eq. ( 4.22] ) would immediately imply a < 0, which is incompatible with 
the known existence of the acceleration period in our recent past (q.e.d.). 

Therefore, the maximum will in general lie ahead in time. After surpassing z = z max 
into the future, the ratio ( 4.14 ) drops fast (in redshift units) to rM(z s ) = — 1, where z s can 



be relatively close to — 1. At this point pu takes the negative value pd{z s ) = ~P°m (l + ; 



with |p£>(£ s )| <C p Q D ■ Prom the previous considerations it is clear that the maximum ( |4.17| ) 
can only exist for 



n% 



>o. 



(4.23) 



u> x (QP x + utf> M )-e(l-n° A ) 

We point out that if the cosmon X would behave as another pure cosmological term (namely 
lux = —1) then e = and the maximum ( 4.17| ) would not exist. On the other hand, for 
v = then again e = 0, and for lux 7^ — 1 the maximum still exists and is located at 



1 + z n 



n A 



u> x n x 



i 



(i/ = 0). 



(4.24) 
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Since lux < this enforces Q\ and Q x to have opposite signs. 

The various scenarios will be further discussed in the next section from a more general 
point of view, where we will link the conditions for the stopping of the Universe's expansion 
with the conditions of existence of the maximum of ( |4.14| ) before reaching the turning point. 
The presence of these two points is a central issue in our discussion of the cosmological 
coincidence problem within the AXCDM . The numerical analysis displaying these features 
is presented in Section 0. But before that, let us take a closer look to the analytical results 
and their interpretation. 

5. Cosmological scenarios in the AXCDM model. Autonomous system. 

For a better understanding of the solution that we have found in Section |3| it will be useful 
to look at the set of equations ( |3.1| ) from the point of view of the qualitative methods for 
systems of differential equations. This should help in the interpretation of the numerical 
analysis performed in Section [?]. As a matter of fact, we will end up by solving quantitatively 
the problem with this alternate method and use it as a cross-check. Again it will simplify 
our analysis if we assume that the cosmon is a barotropic fluid with lux = const. Let us 
introduce the convenient variable 

C = -ln(l + z). (5.1) 

In this variable the remote past (t = 0) lies at £ = — oo whereas the remote future (t = oo) is 
at £ = +oo. In the previous section we have solved the AXCDM cosmological model exactly 
for any value of the spatial curvature, but in practice we will set the curvature parameter Q° K 
to zero to conform with the small value suggested by the present observations M, 0] . This 



will allow us to deal with the basic set of equations (3.1) as an autonomous linear system 
of differential equations in the variable (|5.1|). Indeed, after a straightforward calculation 
we find 

-[ua m + (1 - v)ax] tt x - va m £lA + va m tt c , 



dC 
dHji 



v (a m - ax)^x + va m Vt,\-va m Q, c 



d( 

-rr = («m - oix) ^x + a m Vt,\- a m O c , (5.2) 



where we have defined £l c (z) = p c (z)j fr c . Notice that by virtue of ( 3.13 ) 



n m ( z ) + n A (z) + n x (z) = n c (z). (5.3) 



We have actually used this equation to eleminate Q m from ( |5.2| ) . The matrix of the system 
Q5.2| ) has zero determinant. The eigenvalues read 

Ai = -ax (1 - v) , \ 2 = -a m , A 3 = 0. (5.4) 
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The corresponding eigenvectors (up to a constant) are 



v.ii - ■ ! v I , v 2 






V3 = 1 • ( 55 ) 



If we define the column vector 



"(0 



V i 
/W(0\ 

n A (C) ( 5 - 6 ) 



then the solution of the autonomous system ( |5.2|) reads 

O(C) = Cx vx e Al c + C 2 v 2 e A2 C + C 3 v 3 , (5.7) 

with 

l-n A ^(a m - a x ) ~ ^(a m - ay) Q A - 1/ 

t>l - jz r , 2 - 7- r , O3 - -j • (5.8J 

l-i/ « m - ax (1 - y) a m - ax (1 - z/) 1 - v 



The values of these coefficients are fixed by the boundary condition 

n(c = 0) 



(&x 



n° A I . (5.9) 



One can check that the solution obtained by this method is completely consistent with the 
one obtained in Section for the flat space case. 

Once the solution is expressed in terms of the eigenvalues, as in ( |5.7| ), we can investigate 
the qualitative behaviour of the dynamical system. Let us describe some possible scenarios 
for the phase trajectories. Let us recall that the cosmon barotropic index lies in the interval 



( 4.11 ). Equivalently, 

-5 < a x < 2 < a m (5 > 0) . (5.10) 

The scenarios that we wish to consider are mainly the ones that are compatible with the 
existence of a stopping point in the future Universe expansion. The reason for this is that if 
there is a stopping point then we can argue (see below) that it can provide an explanation 
for the cosmic coincidence problem. Of course, this "stopping point" is actually a "turning 
point" , meaning that the Universe ceases its expansion there, but it subsequently reverses 
its evolution (then contraction) towards the Big Crunch. Among the possible cosmological 
scenarios let us remark the following: 

• I) < ax < 2 and v = 0. In this case the two non-vanishing eigenvalues are both 
negative, Ai < , A2 < 0, meaning that there is an asymptotically stable node. All 
phase trajectories of the system in vector space ( J5lf ) become attracted (as £ —>■ 00) 
towards the vector node 

n}= nU. (5.11) 
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For Q^ > this fixed point is obviously de Sitter space-time. However, if we further 
impose the condition S7^ < 0, the convergence to that node will be stopped because 
the third component of the vector solution (|5.6|), Q c (z) = p c {z)/ p® = H 2 (z)/Hq, 
cannot become negative. The interpretation of this case is clear: we are in a situation 
where we have a quintessence cosmon entity X (—1 < uj x < —1/3), with density (|3.7j), 
superimposed on a strictly constant cosmological term: 

n D (z) = n° A + n° x (i + z) ax . (5.12) 

The requirement Q A < eventually stops the Universe expansion at some point in 
the future because the X density gradually diminishes with the cosmic time and the 
r.h.s. of ( 3.10| ) becomes zero. We remark that the possibility that 09 < should not 



be discarded a priori. Even though it is clearly unfavored by measurements within 
the standard ACDM model g, |, |, §], it is perfectly possible within the AXCDM 
whose generalized cosmic sum rule ( |3.12 ) for flat space reads 



n° M + n° A + n x = i. (5.13) 

For example, if we take the approximate prior Q^ = 0.3 (typically from the LSS 
inventory of matter obtained from galaxy clusters analysis [|]]), then a value Q A < 
can be easily compensated for by a positive cosmon density at present: Q x > 0. 
Notice that the measurements of distant supernovae and CMB should in principle be 
sensitive to the total DE content Qd = $7 a + £l x , and of course to the total matter 
content. Moreover, being a x > we see that Q A < is compatible with the existence 
of the maximum of r(z) at z = z max - cf. Eq.( 4.19| ). 



II) —5 < ax < and v = 0. Here the non-vanishing eigenvalues are of different 



sign: Ai > 0,A2 < 0. It follows that (5.11) becomes a saddle point and the phase 



trajectories depart from it more and more with the evolution. Equation ( p. 12 ) still 
holds in this case, and for Q x > the cosmon density increases indefinitely with the 
expansion. In fact, X behaves here as a standard phantom field (cf. Fig.ffl) and pro- 
duces a super-accelerated expansion beyond the normal de Sitter expansion. In these 
conditions the Universe's evolution strips apart gravitationally bound objects (and 



ultimately all bound systems), i.e. it ends up in a frenzy "Big Rip" [15|. However, if 



£l x < the r.h.s. of ( |3.10 ) eventually becomes zero and the Universe stops its expan- 



sion at some point in the future. We can also see these features in terms of Eq,fl5.7|). 
The sign of the coefficient C\ determines in this case whether there can be stopping 
(Ci < 0) or not (C x > 0). Since u = here, we have C x = 1 - n° A - U° M = n° x 
and therefore it is the sign of Q x that controls the stopping along the lines men- 
tioned above. Notice that although uj x < — 1, the Big Rip is impossible for fl x < 
because the term — p x (1 + 3uix) on the r.h.s. of (2.13) is negative, therefore the 



cosmon in this instance actually collaborates with matter to enhance more and more 
the cosmological deceleration until the Universe reaches a turning point. In this case 
the cosmon X acts as "Phantom Matter" (PM), the only component of the DE that 
still preserves the SEC (see Figl). We stress that even though X behaves as PM 
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• 



here, the sign of the product ax Px is positive and hence the effective behavior of 
the AXCDM model remains quintessence- like - see Eq.( |2.18] ). 

Ill) ax = and v ^ 0. Two eigenvalues are zero (Ai = A3 = 0) and A2 < 0. One can 
check that the following node is the end point of all trajectories: 



n* in (?)= nl-vn% 



(5.14) 



This situation corresponds to effectively having two variable "A's" , one is X (which 
has EOS parameter lox = — 1) and the other is A. Both "A's" are variable because 
v 7^ 0, but the sum pj\{z) + px(z) = const, i.e. Pd{z) behaves as a strictly constant 
cosmological term - cf. Eq. (|2.8|) . With Vfi* M < 1, there is no possibility for stopping 
the expansion under these circumstances, no matter how A varies, because £l c (z) 
(hence H 2 (z)) remains always positive. Although this case is uninteresting for our 
purposes, it shows at least that there is no turning point in the Universe evolution if 
the cosmon has the exact EOS of a cosmological term. In Section H we also proved that 
under these circumstances there is no maximum of r = pn/Pm either. The Universe 
moves asymptotically towards an ever-expanding de Sitter phase characterized by an 
effective cosmological constant given by Q,\ + Q x = 1 — Q^j > 0. We conclude that 
a necessary condition for the existence of a turning point in the AXCDM model is 
that the cosmon EOS parameter is different from —1; in other words, ax = is to 
be excluded from Eq. ( |5.10D , This will be understood implicitly hereafter. 

IV) < ax < 2 and v < 1 (including of course v < 0). As in case I above the two 
eigenvalues are both negative Ai < , A2 < and there is a (^-dependent) node onto 
which all phase trajectories are attracted, namely 

Wv(") = fe • (5-15) 

vfey 

For v = 0, fijy(0) = CIJ and we recover the case I above. Again the attraction 
towards the node can be stopped under suitable conditions. In this case we have to 
require that 

-±- — <0 => Vi\<v<\. (5.16) 

Being ax > we see that this choice is compatible with the existence of the maximum 
at ^max> see Eq.( [4.19| ). In contrast to scenario I, we find that to achieve stopping it 
is not necessary that Q^ < 0. This novel trait appears thanks to the non-vanishing 
interaction between the two components of the DE in the AXCDM, which is gauged 
by the parameter v - as seen from the first two equations in (|3.1| ). Therefore, we 
conclude that for quintessence-like cosmon we can have stopping of the expansion 
not only for Sl\ < but also for < O^ < v < 1. 
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V) < ax < 2 and v = 1. In this case we have the same eigenvalues as in scenario 
III, but the situation is different. The asymptotic node for the trajectories is 

n* v (v) - | n° A - n° M (i - a x /a m ) | . (5.17) 

1 -0^(1 - ax/a m ) 



Using the sum rule ( 3.12 ) it is easy to see that for ax = this fixed point coincides 



with ( 5.14D for v = 1, as it should. In contrast to scenario III, in this case the 



stopping is in principle possible by requiring that the third component of ( |5.17| ) 
becomes negative. This condition yields 

Q° - 1 
ax < -^ a m . (5.18) 

Using a m = 3 for the matter-dominated epoch and the prior O^ = 0.3, we find 
ax < —7 (i.e. u>x < —1-33) to achieve stopping. This would be no problem, were 
it not because the nucleosynthesis constraint derived in Section H forbids that the 
product | fax | is too large as in the present case. So this scenario is actually ruled 
out by nucleosynthesis considerations. 

• VI) < ax < 2 and v > 1. Here one eigenvalue is positive and the other negative 
(Ai > , A2 < 0) as in scenario II. Therefore the point (5.1 1]) is again a saddle 



point from which all trajectories diverge as £ — > 00. This runaway, however, can be 
stopped provided C\ < in ( |5.7[) . Indeed, since the eigenvector vi defines a runaway 
direction the third component of (^J) will eventually become negative, which is the 
condition for stopping. In the absence of the nucleosynthesis constraint ( [4.12 ), the 



and 



condition C\ < would be a bit messy, but thanks to this condition we have in good 
approximation 

Ci = 1 -^_nS £ (^- M ^i- j ji_ 

1 — v uj m — ujx + e 1 — v 

C 2 = ^M^rn-Cx) „ ^ (5 _ 2Q) 

Wm - u x + e 
Thus, within this approximation, the stopping condition for this case is 

d<0 ^ l_ii<O M , (5.21) 

where e - defined in Q4.7| ) - is assumed to be small to satisfy the nucleosynthesis 
constraint. In the present instance, since we assume v > 1 we have to compensate it 
with small enough ax- The condition ( |5.2l| ) is satisfied by choosing 

n° A < 1 + {v - 1) n° M . (5.22) 

This value could still perfectly lie in the preferred range of values for il^ near the 
standard flat ACDM model, i.e. Q A ~ 0.7 |^, ||[ ^, ||. But of course with a dramatic 
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difference: with this same value we can have turning point in the Universe evolution 
within the AXCDM model. Combining the stopping condition ( |5.22| ) with the con- 



dition ( 4.19| ) for the existence of the maximum of r(z) under ax > 0, we see that 



we must have fl A < 1 + (y — 1) £l M < v, hence f^ < 1. Our usual prior for this 
parameter (f^- = 0.3) automatically guarantees it. 

• VII) —5 < ax < and v < 1. In this case we deal with phantom-like cosmon. 
The eigenvalue Ai is positive. Therefore the point (|5.11| ) is a saddle point, similar to 
scenario VI. Again the stopping condition is ( |5 . 2 1| ) . However, here v < 1 and thus 
we have 

n° A > i - (i - v) n° M . (5.23) 

Using the prior < rfi M < 1 we find 

n° A - v > i - n° M - v (i - n° M ) > o. (5.24) 

As ax < the previous equation is compatible with ( |4.19| ), and so the simultaneous 
existence of the stopping and the maximum is guaranteed. 

• VIII) — 5 < ax < and v > 1. Again we stick to phantom- like cosmon here. The 
eigenvalue Ai is nevertheless negative. The trajectories are attracted to a node of the 
same form as ( |5.15|) . However, since v > 1 in this case the stopping condition takes 
the form 

A - - < ^> Vi\>v. (5.25) 

On account of ax < this choice is compatible with the existence of the maximum 
at z max , see Eq.( [4.19| ). 

These are the main scenarios we wanted to stand out. From the previous analysis it is 
clear that de Sitter space-time is not the common attractor for the family of cosmological 
evolutions of the AXCDM model, even in the absence of a turning point. Let us notice the 
following interesting feature: the signs of Qa( z ) and 0,x(z) may change when we evolve 
from the past into the future. This can be easily seen from ( |5.7| ). For example, in scenario 
VI (where v > 1), and requiring stopping (C\ < 0), we have 

Qa(z^oo) ~vC 2 {l + z) am >0, 
Q A (z -» -1) ~ v d (1 + z) ax f 1 "") < , 

n x ( z -► oo) ~ - uC2fXm (i + z) am < o , 

«m - ax 
n x {z -» -1) s (1 - v) d (1 + z) ax {1 - v) > , (5.26) 



with C2 — Q° M > 0. We have also used a m > ax - see Eq. ( 5.10| ), The sign of the 



cosmological term goes from positive in the past to negative in the future, and this is the 
physical reason why the Universe halts in this case. On the contrary the cosmon density 
in this same scenario is seen to be negative in the past and positive in the future. At the 
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turning point z s , where the expansion rate H(z) vanishes, the sum Qd(z) = Qa(z) + Qx(z) 
takes on the negative value 

n D ( Zs ) = -n° M (i + z s f , (5.27) 

and therefore it cancels against the matter contribution 7 . In contrast, in scenario VII 
(specifically for the situation v < 0) the sign changes of Q\(z) and £lx(z) from the past 
to the future are just opposite to scenario VI, as can be easily checked from the formulae 
above, so that here the cause for the stopping is actually attributed to the negative energy 
of the cosmon in the future (rather than to the effect of a negative cosmological term). 
These are two complementary situations. In both cases the cosmon field (or the effective 
entity it represents) has negative energy either in the past or in the future. It thus behaves 
in these cases as a phantom-like field with negative energy, i.e. as Phantom Matter (cf. 
Fig. [I]). However, as in scenario II above, the effective behavior of the AXCDM model 
remains quintessence-like because ax Px > 0. 

6. Crossing the u e = — 1 divide? 



The recent analyses [13, 14] of the observational data imply the possibility of an interesting 
phenomenon known as the crossing of the cosmological constant boundary fl5f . The essence 
of this phenomenon is that the parameter of the dark energy equation of state, uj e = uj e (z) 
at some positive redshift z = z* attains, and even trespasses, the value u> e (z*) = —1, i.e. 
the function u e (z) crosses the cosmological constant divide uj\ = — 1. Some analyses of the 
cosmological data mildly favor the transition from the quintessence-like to the phantom- 
like regime with the expansion at a small positive redshift of order 1. Let us recall from the 
analysis of Ref. |37], |38| that in models with variable cosmological parameters we generally 
expect an effective crossing of the cosmological constant barrier at some point near our 
time (whether in the recent past or in the near future). This crossing actually refers to the 
effective EOS parameter u e = uj e (z) associated to the model description in the effective 



"DE picture" [38], where the DE is self-conserved. This picture is in contradistinction 



to the original cosmological picture where A is either exchanging energy with matter or 



is varying in combination with Newton's gravitational constant G[38[. Of course in the 
original picture oj\ = — 1, whether A is strictly constant or variable. The crossing, therefore, 
is visible only in the DE picture, which is the one assumed in the common parametrizations 
of the DE, and therefore the one used in the analysis and interpretation of data. 

What is the situation with the AXCDM model? In the present instance we do not 
necessarily expect such crossing. The reason is that the AXCDM model is not a pure 
A model (variable or not) due to the presence of the cosmon; in fact, A and X form 
a non-barotropic DE fluid with uj e ^ — 1 from the very beginning. For this reason the 



AXCDM model does not necessarily follow the consequences of the theorem [38]. It may 



7 As warned in Sect, ti the formal structure of ( pi.lq ) should not lead us to think that pn — Pd(z) remains 
necessarily positive in this model for the entire history of the Universe. The effective EOS parameter 



Lj e = u> s (z) of the AXCDM model (2.7) is actually singular at a point zq > z s (i.e. before stopping) where 



Pd(zo) = 0. Of course this is a spurious singularity because all density functions remain finite at z = zq. 



After that point (i.e. for z s ^ z < zo) we have £Id{z) < 0, and at z = z s we have (5.27) exactly. 
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exhibit crossing of cv e = — 1, but it is not absolutely necessary. As a matter of fact, in 
certain regions of the parameter space there is no such crossing at all, but in other regions 
there is indeed crossing. In the following we concentrate again only on the flat space 
case. We shall prove next that whenever there is crossing of uo e = —1 within the AXCDM 
model it can only be of a specific type in order to be compatible with the existence of a 
turning point. Specifically, the crossing should go from phantom-like behavior in the past 
into quintessence-like behavior at present or in the future. It can never be the other way 
around (if a turning point in the Universe evolution must exist). 



The proof goes as follows. From Eq. (2/7) we have (1 + u e )pD = (1 + ^x) Px- The 



explicit solutions for the individual DE densities Pd(z) and px{z) found in Sections |3| and 
|5| - see e.g. Eq. ( |3.8| )- reveal that these functions are nonsingular for finite values of the 
scale factor, i.e. for the redshift values in the entire range — 1 < z < oo. Therefore, in 
order to have uj e (z) = — 1 at the transition point z = z*, the cosmon density must vanish 
at that point: px(z*) = 0. The explicit expression for ilx(z) = Px{z)/p° c ( c f- Eq. ( |3.8| )) 
can be rearranged in the form 



^1 = (i + 6)(i + zfx (i-") - 6(1 + *)«™ , 

A 



b(l + z) 



+ n + z) ax i 1 -")-^ 



where 



i/(i + w ro ) n° m u(i + u m ) n° m 



(6.1) 



(6.2) 



1 + u m - (1 + wx)(l -v) Sl Q x Lu m - uj x + e Q° x ' 

For v = it is patent that the expression (pll) shrinks to fl3.7| ) and it does never vanish 
for z > —1. In contrast, for v ^ the two terms in ( |6.l| ) may conspire to produce a zero. 
From the structure of Qx( z ) it is clear that 0,x(z) can have at most one zero in the entire 
interval — 1 < z < oo, namely at z = z* given by 



] a X (l-v)-a m 

1 + ,- = { — } . (6.3) 

Consequently, the crossing of the io e = — 1 line can happen at most at one value in this 
redshift interval. Let us now concentrate on the case when the DE evolves with the ex- 
pansion of the universe from the quintessence regime into the phantom regime at some 
positive redshift z* . In this case for z > z* we have cu e > — 1 whereas for z < z* we have 
uj e < — 1 and the present value of co e is phantom-like, i.e. uJ e (0) < — 1. Since by assumption 
the crossing happened in the past at some z* > 0, there can be no crossing in the future, 
since there is at most one crossing. It follows that in the future io e remains below — 1, and 
since the present DE density is positive in our model (p° D > 0), the function p£> = pd(z) 
can only grow to larger and larger positive values with the expansion of the Universe - see 
Eq. ( [2,16 ). Therefore the expansion rate H given by ( [2.12 ) will also keep growing forever, 



and hence no stopping is possible. This is tantamount to say that uj e will exhibit no singu- 
larities in the future, which is a necessary condition to trigger stopping of the expansion. 
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We conclude that this type of cosmological constant boundary crossing is not compatible 
with the condition of stopping of the expansion (q.e.d.) 

On the other hand the alternate transition from the phantom-like to the quintessence- 
like regime with the expansion of the Universe is characterized by u e (0) > —1, for z* > 
0. This type of crossing in the past (z* > 0) is perfectly compatible with the stopping 
condition, and we can ascertain analytically under what conditions it takes place. From 
Eq. ([T^) we see that in principle the crossing can only take place under one of the two 
following conditions: 

(i) b>0 and ax (1 — v) — a m < 0, or (ii) b < — 1 and ax (1 — v) — a m > . (6.4) 

However, option (ii) is impossible because the restriction (|4.11| ) on the range of u>x and 



the nucleosynthesis condition ( |4.12[) imply together that 



a x (1 - v) - a m = 3 (w x - u m - e) < . (6.5) 

Hence we are left with the condition b > 0. From ([T^) we finally conclude that the desired 
crossing in the past will take place whenever Q x and v have the same sign: 

(n° x > 0, v > 0) or (n x < 0, v < 0) . (6.6) 



Furthermore, on comparing the condition for the existence of the crossing ( |6.3| ) with the 
stopping conditions listed in Section |5| we see that they are compatible. Explicit examples 
of this kind of crossing are illustrated in the numerical analysis of Section 0. 

As stated above, some analyses of the cosmological data mildly favor the transition 
from the quintessence-like to the phantom-like regime, which is opposite to the kind of 
behavior we find for the present model (if the existence of both a transition point and a 
turning point is demanded). However, the crossing phenomenon has never been established 
firmly to date, and the more recent analyses [jl | cast some doubts on its existence. Actually, 
the same kind of analyses leave open the possibility that the crossing feature, if it is there 
at all, can be of both types provided it is sufficiently mild; that is to say, we can not exclude 
that in the recent past the Universe moved gently (i.e. almost tangent to the to e = — 1 
line) from an effective phantom regime into an effective quintessence regime. In such 
case the average value of oj e {z) around z = z* would stay essentially u e = — 1 (therefore 
mimicking the standard ACDM). Interestingly enough, the AXCDM model can perfectly 
accommodate this possibility (see next section for explicit numerical examples), but at the 
same time the model is also perfectly comfortable with the absence of any crossing at all 
while still maintaining a turning point in its evolution. 

7. Numerical analysis of the AXCDM model 

7.1 Parameter space 

In the previous sections we have solved analytically the AXCDM cosmological model and 
we have performed a detailed analysis of its analytical and qualitative properties. In 
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the present section we want to exemplify these properties by considering some numerical 
examples. In the first place it is important to check that the region of parameter space 
that fulfills all the necessary conditions is a significant part of the full parameter space of 
the AXCDM model. The latter is in principle defined by the following 5 parameters 

(n° M ,n° A ,n° x ,u; x ,u). (7.1) 

However, we will take a prior on £79, because the matter content of the Universe can be 
accounted for from (galaxy) clusters data alone Q, i.e. from LSS observations independent 
of the value of the DE measured from, say, supernovae or CMB. The approximate current 
value is |j|]: f^ = 0.3. Moreover, since we want to restrict ourselves to the flat space case 
(Qjs: = 0), we find from the sum rule ( 3.1 2| ) that 



n° D = n° A + n° x = o.7. (7.2) 

This choice of the cosmological parameters will be implicit in all our numerical analyses. 
We can then take e.g. S7^ as an independent parameter. Thus in practice we are led to 
deal with the basic three-dimensional parameter space 

(n° A ,u x ,u). (7.3) 

It is on this smaller parameter space that we must impose the necessary conditions to 
define the "physical region" of our interest. These conditions are essentially three: 

• i) The nucleosynthesis bound (cf Section ||) 

M<10%, (7.4) 

i.e. the condition that the DE density is sufficiently small as compared to the matter 
density at the nucleosynthesis time; 

• ii) The stopping of the expansion, i.e. the condition of reaching a stopping point 
in the Universe evolution from which it can subsequently reverse its motion (see 
Section ||); 



• 



iii) The condition that the ratio of the DE energy density to the matter density - see 
Eq. ( [4 .![ ) - is bounded and stays relatively small throughout the entire history of the 
Universe, say 

— < 10 (-KK+00), (7.5) 

tq being its present value (see Eq. ( f4.4[ ) and Section || in general). 



The "physical region" is defined to be the subregion of the parameter space ( |7.3[ ) 
satisfying these three conditions at the same time. It is illustrated in the 3-D-plot in Fig. [2]. 
In Fig. |3| we can see the projection of the full 3D physical region onto three orthogonal 
planes Q\ = const., ujx = const, and v = const, situated outside this region. It is also 
convenient to show cross-sections of the physical region, i.e. slices of the relevant volume 
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Figure 2: The 3D-dimensional "physical region" of the AXCDM model, consisting of the subset 
of points ( |T.3| ) satisfying the three requirements: i) Nucleosynthesis bound (7.4), ii) Existence of a 
turning point in the Universe evolution, and iii) Ratio of DE density to matter density, r — Pd/ pm. 



under the bound (7.5) in the entire Universe lifetime . 



for fixed values of the other parameters, just to check that the physical region is not just 



formed by a closed 2D-surface, but rather by a 3-D-volume contained in the space (7.3). In 
Fig.|| we show three slices Q\ = const., u>x = const, and v = const, for particular values 
of these parameters that illustrate this feature. We note that in Fig.[3| we have allowed for 



uix to be slightly out in the range (4.11), only for completeness. In fact the region that 
does not satisfy condition ( 4.11| ) is seen to become narrower the larger becomes ujx, and 
therefore it is rather restricted. 



7.2 Expansion rate and effective EOS 

In Fig.|5] we show the expansion rate H = H(z) of the AXCDM model as a function of the 
redshift for given values of the parameters: lux = —1.85, il A = +0.75, v = —vo, where vq 
is defined in ( 2.22J ). This selection of parameters corresponds to scenario VII in Section 
|5l Since Eq. ( 5.23^ is fulfilled, there must be a stopping point ahead in time, which is 
indeed confirmed in Fig. [5] where we see that H vanishes at some point in the future. As 
discussed in Section || the cause for the stopping in this case is actually attributed to the 
negative energy of the cosmon in the future (i.e. to the presence of "Phantom Matter" 
rather than to the effect of a negative cosmological term). The presence of a stopping point 
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Figure 3: The projections of the 3-D-dimensional physical region of Fig.pl onto the three planes: 
(a) n\ = const., (b) u>x = const, and (c) v = const. 



implies the existence of a maximum for the ratio between DE and matter densities, and this 
fact provides a natural explanation for the coincidence problem. In the present case the 
maximum exists because Eq. ( |5.24| ) is satisfied. Notice that the choice of parameters used in 
Fig. P has the attractive property that the value of Q\ is very close to the current one in the 
standard ACDM model. However, it should be remembered that in all numerical examples 
discussed in this section we have fixed O^, as in (|7.2| ). This means that in the particular 
numerical example illustrated in Fig. |5J, the value of the cosmon density at present is rather 
small and negative: Q^ = —0.05. 

Very important for the future program of cosmological experiments J?], ||, [| is the 
issue of the equation of state (EOS) analysis of the dark energy from the observational 
data. In the standard description of the DE, it is assumed that the DE consists of a 
self-conserved entity, or various entities of the same kind which produce a conserved total 
DE. As a particular case we have the standard ACDM model, where the DE is just the 
constant A. But we can also have simple quintessence models, with a single scalar field or 



multicomponent scalar fields [19, 27 1. In the latter case one of the fields can be a (standard) 
phantom (cf. Fig.[l|) and this has been proposed as a mechanism to explain the crossing of 
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Figure 4: Different slices of the 3D volume representing the physical region, see Fig.pl. The clotted 
areas fulfill the nucleosynthesis condition. The shaded regions satisfy both the nucleosynthesis and 
stopping conditions. If we further require the condition (7.1) on the ratio r = po/pmi we are 
left with the darkest shaded region. Slices shown in the figure: (a) Slice f2° = 0.75; (b) Slice 
lox = —1-85; (c) Slice v = —Vq, where vq is defined in ( 2.22 ). 



the cosmological constant divide [15]. In practice the experimental EOS 



PD = w exp Pd 



(7.6) 



has to be determined by observations. This equation will be confronted with the theoretical 
predictions of any given model. For the theoretical analysis it is convenient to consider the 
general notion of effective EOS of the DE |3(|], namely the EOS that describes dark energy 
as if it would be a single self-conserved entity without interaction with matter . We write 
it as 

p D = uj e (z)p D , (7.7) 

where cu e = ui e (z) is the theoretical EOS parameter or function in the effective DE picture, 



which is the one that should be compared with ui exp in ( |7.6|) , In the original model the DE 
could actually be a mixture of entities of very different kind, as indicated in ([2,1|). This is 
indeed the case for the AXCDM model where we have a non-barotropic mixture of cosmon 
and running cosmological constant. From the explicit solution of the model obtained in 
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Figure 5: The Hubble function H = H(z) for the values of the parameters: lux 
0.75, j/ = —i/o. This situation illustrates an example of Scenario VII (see Section H). Condition 
( 5.23Q is satisfied and hence there is a stopping point in the future where H vanishes and the 
expansion of the Universe stops. 



Section |3| we can study the numerical behavior of the effective EOS of the AXCDM model, 
see Eq. Q . 

Before displaying our EOS results, we recall that the most recent WMAP data, in 
combination with large-scale structure and supernovae data, lead to the following interval 
of values for the EOS parameter ||: 



u, 



cxp 



-1.06 



-0.13 
-0.08 



(7.8) 



This result is quite stringent, and it does not depend on the assumption that the Universe 
is flat. However, it does depend on the assumption that the EOS parameter does not evolve 
with time or redshift. In this sense it is not directly applicable to the effective EOS of the 
AXCDM model. In general the experimental EOS parameter u; e xp need not be a constant, 
and one expects it to be a function of the redshift, uj exp = u e ^ v {z). Since this function is 
unknown, one usually assumes that it can be parametrized as a polynomial of z of given 
order. This is, however, inconvenient to fit cosmological data at high redshifts; say, data 
from the CMB. And for this reason other parametrizations have been proposed and tested 
in the literature [11, O, 13]. As an example, consider the following one [O], 



0J e (z) = OJQ +LOq 



(1+*) 



0=1,2) 



(7.9) 



where u> = (dco e /dz) z= o. The last parameter contemplates the possibility of a residual 
evolution of the DE even if cuq = —1. Notice that the asymptotic behavior of co e (z) 
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Figure 6: The effective EOS parameter of the AXCDM model, uj e -defined in Eq. (2/7) -, as a 
function of the cosmological redshift, for different scenarios: (a) u>x = —1-85, £1° = 0.75, v = —Vq] 
(b) As in (a) but for v = 0; (c) As in (a), but for v = +v ; (d) ux = -0.93, ft^ = -2, v = +0.96. 



in this parametrization is tamed: for p = 1 one has u; e (+oo) = ujq + lo' , and for p = 2, 
w e (+oo) = cjo- This kind of parametrizations can be useful for a simplified treatment of the 
data. However, in general one should not be misled by them. The underlying dark energy 
EOS can be more complicated and may not adapt at all to a given parametrization. In this 
respect it should be kept in mind that there is also the possibility that the experimental 
EOS (|7.6|) is mimicked by other forms of DE that look as dynamical scalar fields without 



being really so. This possibility has been demonstrated in concrete examples in [37] and 
further generalized in 



]. In particular, in the analysis of [37] it is shown that the "effective 
EOS of a running cosmological model" can be analytically rather complicated; it does not 
adapt to the form ([O]) and nevertheless it provides a good qualitative description of the 
experimentally fitted EOS from the data [O]. In the following we will show, with specific 
numerical examples, that the effective EOS of the AXCDM model may account for the 
observed behavior of the DE in the redshift range relevant to supernovae experiments. In 



Section 7.4 we will show that even the asymptotic behavior of the AXCDM model at very 



high redshifts (relevant to the CMB analysis) completely eludes the class of models covered 



by simple parametrizations like (7^9). These parametrizations are, therefore, not model- 
independent, and the fits to the dark energy EOS derived from them should be handled 
with care and due limitation. 
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Figure 7: Detail of the discontinuity of the EOS function u) e = u) e {z) at a point z$ < —0.638 in 
the future, corresponding to case (a) of Fig. g. At this point pd(zo) = 0. 



In Fig.[6| we plot o> e = 0J e (z) for the AXCDM model as a function of z for some 
particular choices of the basic set of parameters ( [7.3| ) and within the typical accessible 
supernovae range j?]]. On comparing these results with the experimental fit (7.8) we should 
stress again that the latter is only valid for a strictly constant EOS parameter. Even so, 
scenarios (a)-(c) of Fig. |6] satisfy this bound already at la. Scenario (d) satisfies it only 
at ~ 2.5a, but here uj e varies substantially with z. 

All the plots in Fig. |(| exhibit a discontinuity at a certain point zq in the future, before 
reaching the stopping point z s . A magnified view of this discontinuity for the case of Fig. 
|6|a can be appreciated in Fig. |7[ The discontinuity occurs when uj e changes sign at a 
point zq where pd{zq) = 0. At this point uJ e {zo) becomes undefined, see Eq. (p7?|). In a 
neighborhood z = zq ±e of the point zq we have co e (z ) = ±oo for e — ► 0. The singularity 
is not a real one, because all density functions remain finite at that point, so its existence 
only reminds us that the description in terms of the EOS parameter is inadequate near 
that point. 

Looking at Fig. || we can see diverse scenarios. For example, in Fig. |6|a we have a 
concrete example of the scenario VII (with v < 0) defined in Section ||, the same one used 
in Fig. ||, namely ujx = — 1.85, fl^ = 0.75, v = —vq. It corresponds to have 0° very near 
the standard ACDM value with the v parameter fixed at the canonical value ( 2.22| ). At 
the same time the cosmon field, X, is phantom-like and it behaves as Phantom Matter in 
the future. We can also appreciate a mild transition of ui e from phantom to quintessence 
behaviour at some point z* near past. This is the kind of possible transition we have 
foreseen analytically for the AXCDM model in Section ||. We can check that the crossing 
condition (i) in Eq. fl6,4| ) is satisfied. If the transition is sufficiently mild and the average of 
uj e around the transition point z* is close to —1 (as it is indeed the case here) it cannot be 
excluded by the present data. Notice that the phantom regime starts to be significant only 
for z > 1.5 in this case. At the moment we do not have enough statistics on this region, 
and it is not obvious that we will ever have, at least using supernovae data alone. Recall 
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Figure 8: The total and individual DE densities, flo = ^d(z), Qa — ^a(z) and fix — £lx{z) 
for two different situations where the cosmon barotropic index is of phantom- type (u>x < — 1) 
and of quintessence-type (lux > —1) respectively: (a) ux = —1.85, D, A — 0.75, v = —Vq) (b) 
u>X = —0.93, fl A = —2, v = 0.96. The two examples correspond to cases (a) and (d) of Fig.g. 

that SNAP 0| is scheduled to possibly reach up to z = 1.7, with a very meager statistics 
foreseen around this high redhshift end. In this sense, a mild phantom behavior around 
that distant regions will be difficult to exclude; in particular, a kind of smooth transition 
of the type Phantom — ► Quintessence (for increasing z) will always be a possibility to keep 
in mind. Notice that in all cases shown in Fig. S the value of u e (z) at z = is very close 
to — 1, namely cu e { z ) = —0.94 in the first three cases, and oo e { z ) = —0.73 in the last one. 
As we have said, these values are sufficiently close to —1 to be acceptable in the light of 
the recent data and the priors used to fit it. 

In Fig. yb we display a situation corresponding to v = 0. This is an example of scenario 
II in Section ||. In this case the cosmological term A remains strictly constant and the 
stopping of the Universe is due to the super-deceleration produced by the cosmon (mainly 
in the future time stretch). Again O^ < 0, but here the cosmon behaves as Phantom 
Matter for the entire span of the Universe lifetime. For v = 0, we can see that in the past 
uj e approaches fast to —1 asymptotically. The reason for this can be seen e.g. in Eq. ( 5. 26] ) 
- and more generally from (6.1) - which shows that in the asymptotic past the cosmon 
density £lx(z) is proportional to v. Figures |6|c,d display other interesting situations. They 
constitute examples of scenarios VII and IV respectively (both with v > 0). In particular, 
in Fig. gc the EOS parameter approaches u> e = — 1 for a long stretch of the accessible 
redshift range by SNAP, and therefore it mimics a constant cosmological term. This case 
is different from the one in Fig.|6]b, where A is strictly constant but the DE is not; here 
the AXCDM model behaves as the standard ACDM model for sufficiently large redshift. 
Finally, in ||d we have a situation where the cosmon is quintessence- like (uj > — 1) and 
nevertheless the effective behavior of the total DE displays a transition from phantom to 
quintessence regime. Again we can check that the crossing condition (i) in Eq. ( |6.4| ) is 
satisfied. 



7.3 Evolution of the DE density and its components 

The behavior of the effective EOS parameter to e is related to that of the DE densities in 
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Eq. ([2,7l). The interaction of the cosmon and the cosmological term produces a curious 
effect; namely, the evolution of the cosmon density, Ox, does not follow the expectations 
that we would usually have for a fluid X carrying a value uj > — 1 or lox < — 1 of its 
barotropic index. For a better understanding of the situation we have plotted in Fig. ||a,b 
the total density Qd(z), as well as the two individual components £l\(z) and Ox (z), for 
two particular cases: one with the same parameters as in Fig. |B|a (in which the cosmon has 
cox < —1 as if it were a standard phantom field), and the other for the same parameters 
as in Fig. gd (where the cosmon has u>x > —1 as if it were a conventional quintessence 
field). We can appreciate from these figures the sign changes of the various DE densities 
with the Universe evolution. In particular, in Fig. ga the cosmon density appears to be 
Ox > in the intermediate past but it starts behaving as Phantom Matter (Ox < 0, cf. 
Fig. ||) near our recent past and also at present, and this trend is enhanced in the future. In 
contrast, for the case of Fig. ||b the X density changes fast from Ox < to Ox > at some 
point in the past (and stays so all the way into the future). Now, due to its interaction 
with the variable A, the quintessence- type cosmon (ux ^ —1) of Fig.||b does not decay 
with the expansion; it actually grows fast as if it were a true phantom! The overall result 
is that, thanks to the concomitant depletion of the A density into the negative energy 
region (very evident in the figure) , there is a net deceleration of the expansion and a final 
arrival at the stopping point. Notice that, in contrast to the hectic evolution of the two 
component Oa(z) and fijf(z), the total DE density Od(^) evolves quite slowly. A different 
behavior is observed for the phantom- type cosmon (u>x < — 1) i n Fig.|8]a. In this case, 
rather than increasing rapidly with the expansion, the cosmon density actually decreases 
slowly (it actually behaves for a while as standard quintessence) and then at some future 
point near our time it "transmutes" into Phantom Matter and causes a fast deceleration of 
the expansion, with an eventual stopping of the Universe. In this case, the total DE density 
Ox)(z) again evolves very mildly. The behaviors in Fig.|||a,b are an immediate consequence 
of Eq. ( 2.18| ), which implies that the AXCDM model effectively behaves as quintessence if 



and only if ax o,nd Ox have the same signs, as it is indeed the case here. 

Recall that in all cases considered in our numerical analysis the total DE density is 
normalized to (|7.2| ) at z = 0, as can be checked in the examples presented in Fig.|8|. In 
this figure it is remarkable the stability of Q,£,{z) in the entire range relevant to SNAP 
experiments [0j both for small v (Fig. ga) and large v (Fig. ||b). This is due to the nucleo- 
synthesis constraint and to the small values of the cosmon parameters ax and/or Q x in 
these examples. To see this, let us first expand the total DE density ( 3.11[) in the small 
parameter e in first order: 

l- V LO m -LO X \ 1-V 



eO° 
+ € m -36 

u m — Ux 






ln(l + zU (l + z) ax . (7.10) 



Next, expanding this equation linearly in z, the ^-dependence cancels and we find a very 
simple result: 

D (z)~0?) + axO x z. (7.11) 
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This expression coincides with the expansion, up to the linear term in z, of the function 



(5.12). The latter is the DE density for the v = case - see scenarios I and II of Sec- 
tion^. Therefore, we arrive at the following remarkable result: due to the nucleosynthesis 
constraint, the DE density of the AXCDM model behaves always (i.e. for any v) as in the 
v = case in the region of small z < 1 (up to terms of order e 2 ). As a consequence the 
evolution of Qd with the redshift will be mild in this region provided ax is small or Q, x 
is small, or both. Since v has dropped in the linear approximation ( |7. 1 1| ) , this result is 
independent of v. Therefore, v can be large (meaning \v\ ~ 1) and still have a sufficiently 
"flat" evolution of Qd( z ) for z < l 8 . We can appreciate particularly well this feature in 
the situation depicted in Fig. ||a, where the cosmon energy density is very low (in absolute 
value) both in the intermediate past and also near the present. Thus both the total DE 
density SId{z) and the cosmological term density, £l\(z), remain virtually constant for a 
long period in our past and in our future - see the plateau tip — const, in Fig. ||a. Although 
in this example the cosmon density, £lx, remains very small in our past and around our 
present ( \£lx\ "C ^a), in the future Qx dives deeply into the Phantom Matter region and 
eventually controls the fate of the Universe's evolution. So the cosmon density, which is 
negligible for the redshift interval relevant to SNAP measurements, does become important 
in the future and is finally responsible for the stopping of the Universe. Similarly, in Fig.|8|b 
we have a large value of v around 1 and Q,n stays also rather flat. In this case the flat 
behavior is due to ax < 1. In both situations represented in Fig.g] the AXCDM model 
mimics to a large extent the standard ACDM model. 

Recall from Section || that in the limit of vanishingly small v and Q. x the AXCDM 
model boils down exactly to the standard ACDM model. The result proven above, however, 
tells us something new: namely, if z is not too large (z < 1, still within most of the SNAP 
range), the coincidence between the two models can still be maintained under appropriate 
circumstances (small ax and/or £l x ), even if v ^ 0, by virtue of the nucleosynthesis con- 
straint on the parameter e. Therefore, the following question is now in order. Under these 
circumstances, how can one possibly differentiate between the two models in a practical 
way (e.g. using the future SNAP experiments)? The answer is straightforwardly simple: 
by looking closely at the EOS behavior which, contrary to the DE density, can be very 
sensitive to the value of v\ In fact, the size and sign of v can play a key role to modulate the 
quintessence/phantom behavior effectively developed by the AXCDM model at moderate 
values of z accessible to experiment. The EOS behavior has been displayed in Fig.|(| but 
the z/-sensitivity of this behavior can be much better appreciated in the comparative set of 
EOS curves shown in Fig. ^, which have been performed for various values of v (with both 
signs) and for fixed values of the other parameters. For a cosmon of the "wx < — 1 type", 
we can see that a positive value of v helps to lift the EOS curve into the quintessence region 
and, therefore, it helps to resemble more to some recent EOS analyses of the data, which 



seem to favor the existence of a mild quintessence-regime in our recent past [11, O]. For 



v < (and small), instead, the curve gently crosses the cosmological constant divide and 



8 This situation is in contradistinction to the RG model presented in Ref. J23J], where the departure of 
the DE density (pa density in that case) with respect to its constant value in the standard ACDM model 
depends only on v. In that model, we must have \v\ <§C 1 if the DE evolution is to be mild. 
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Figure 9: Comparison of the effective EOS parameter of the AXCDM model, uj e , for fixed values 
of lox and £1^ (i^x = —1.85, fi° = 0.75, as in Fig. 0a) and different values of v in units of vo, 
Eq. (2.22). All curves give w e (0) = —0.94 at the present time. 



plunges slowly into the phantom regime at higher and higher z. For a "uj > — 1 type" 
cosmon, instead, these features are opposite (cf. Fig.^Jd). Finally, for v = the curve 
stabilizes very rapidly (at increasingly high redshifts) along the flat line u) e = — 1 and the 
AXCDM model looks in this case essentially as the standard ACDM both from the point of 
view of EOS behavior and constancy of the DE density. This is because ax < in Fig.|9|, 
meaning that Qd( z ) -^ ^A f° r z sufficiently deep in the past (see next section). But even 
in this case the effective EOS starts changing softly at some point near our past, and then 
steadily fast in the future, until it breaks into a singular (but spurious) behavior at a point 
zq where pd(zq) = and u> e (zo) becomes undefined. 



7.4 Asymptotic regime of the EOS in the past 

Let us study the behavior of the effective EOS in the asymptotic past, z S> 1. This will 
help to substantiate the behavior of the numerical examples considered above and reveal 
additional features. From Eq. (|6.1[) we find 



n x (z > 1) = < 



-bn x {l + z) a ™, fori/ ^0 



for v = 0, ax < 



(7.12) 



n° x (1 + z) ax for v = 0, a x > . 
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where b is given by (|6.2| ). Consider next the exact expression for Cljj(z), given in ( 3,11| ). 
The asymptotic behavior of this expression in the past is 



Q D (z » 1) = < 



U)m— (^X+e 



n° m (i + zy 



n% 



for v ^ 

for v = 0, ax < 



(7.13) 



Hi (1 + z) ax 



l x yx -r <*) - for i/ = 0, ax > . 

For completeness let us also quote 



fi A (z » 1) = < 



nl 



(7.14) 



for z/ = . 



We can indeed check out the predicted trends for the various DE densities on inspection 
of Fig. p. Substituting equations ( 7.12 ) and ( |7.13| ) into fl2.7| ), and using ( |6.2j ), we find the 
corresponding asymptotic behavior of the effective EOS parameter in the past: 



w e (z> 1) =-! + {! + UJ X ) 



tt D (z > 1) 



oj m for v 7^ 
— 1 for i/ = 0, ax < 
k ujx for zv = 0, ax > . 



(7.15) 



On the one hand the asymptotic results for u> e in the case v = were easily foreseeable 
from the structure of Eq. ( |2.7| ), because in this case (that is, in the absence of interaction 
between the cosmon and A) the growing and decreasing properties of the cosmon density 
strictly follow its canonical quintessence or phantom character in each case. In other words, 
for !/ = fla cosmon carrying a phantom-type or quintessence-type wx-label does properly 
behave as a standard phantom or quintessence, respectively. (Not so, however, when v ^ 0, 
as we have seen!) On the other hand the asymptotic result for u e in the case v ^ came 
a bit of a surprise. It implies that at very high redshift the effective EOS of the DE in the 
AXCDM model coincides with that of matter-radiation. This does not cause any harm to 
nucleosynthesis, for as shown by Eq. ( 7.13j ) the asymptotic DE density remains bounded by 
the constraint ( 4.10Q . This prediction of the AXCDM opens the possibility that the CMB 
data could actually be blind to the EOS analysis of the DE. Indeed, as we have mentioned 
above the CMB fits usually assume parametrizations of the form ( |7.9[ ) and the like. In 
general one replaces Eq. ( 2.16 ) into the expansion rate ( |3.10| ) to get (for £l° K = 0) 



H 2 (z 



m 



n° m (1 + z) am + n° D exp <^ 3 / dz 



f l + u e (z') 



l + z> 



(7.16) 



Whatever it be the "well-behaved" parametrization assumed for to e at high redshift, if it 
leads to powers of 1 + z in the function H = H(z), the fact that z — 1100 for the CMB 
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analysis will unavoidably provide a fit of to e very close to 
Eq. ( fT!J| ) for p = 1, we find 



■1. For example, in the case of 



H\z) = Hl \S&(l + z) a 



+ n° 



D 



[l + Z 



,3(l+o;o+Wo) 



exp 



"3Wn 



l + z 



(7.17) 



This equation reduces to the standard ACDM one for uq = — l> w o = ^> as expected. At 



the same time it is clear that fitting data at very high redshift (e.g. CMB data) with ( 7.17| ) 
will naturally convey uj e ~ — 1 - if one wants to match at the same time the LSS and 
supernovae data both suggesting that Q° D is of the same order as &° M - Therefore, finding 
u> e ~ — 1 from (7.17) can be an intrinsic bias or artifact of the parametrization itself. 



In contrast, if the underlying theory of the DE would exhibit a behavior of the AXCDM 
type ( [7.15 ), the CMB data could not be sensitive to the EOS parameter in its conventional 
form and the observations should be interpreted in a different way. Namely, the canonical 
result <jj e ~ — 1 could only be derived for "small" or intermediate redshift of the order 
of the distant Type la supernovae, where the AXCDM model can resemble to a large 



extent the standard ACDM. In the new scenario the result ( 7.15 ) for u ^ could still be 
detected by observing a "mass renormalization" effect of the cosmological parameter fi^ 
when comparing the two sets of results for small and high redshifts, i.e. when comparing 
supernovae data fits with CMB fits. Indeed, substituting ( 7.13J) in ( |3.10| ) one gets, at very 
high redshift (neglecting the curvature term) 9 : 



H 2 



H 2 n° m (i + z) 



where the "renormalized" cosmological mass parameter is: 



n° 



o 



u>„ 



cox +e 



(7.18) 



(7.19) 



The relative difference |f^ — f2^J/f2^ between the two determinations of O^ is just given 
by the nucleosynthesis constraint on the parameter r = po/pm, see Eq. ( [4.10 ) . Thus it 
could lead to a measurable ~ 10% effect that would provide a distinctive signature of the 
model. There is nothing odd in the DE behaving this way; the DE density ( 7.13J ) is still 
there all the time even though the conventional EOS analysis could miss it. In this case 
one must resort to a different interpretation which may lead to its determination. 

Let us also note that since uj m ^ 0, then in those {y ^ 0) cases where there is a crossing 
point z* from phantom into quintessence regime there must necessarily be a point in the 
past where pr> becomes zero and changes sign. This produces another spurious singularity 
in the function oj e . (Recall from Section |B] that we cannot have more than one crossing 
point z* .) In this case the evolution of uj e from the asymptotic regime u) e {z 3> 1) ~ u m > 
up to the singularity, is such that u> e stays above — 1. After the fake singularity is crossed 
towards z = z*, the value of io e emerging from it is smaller than —1, and remains so until 
reaching the crossing point z* where again u e becomes larger than —1. 



9 We note from the exact formula (3.11) that the total DE density in the AXCDM does not scale in 
general as matter-radiation ~ (l + z) Qm , because of the additional term ~ (1 + z) s ^ 1+u ' x ~' 1 ' which plays an 
important role for z ~ 1. It is only at very high redshift that the aforementioned scaling takes place. 
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Figure 10: The ratio r = pd/ Pm, Eq. (1-1), as a function of the redshift (in units of the current 
ratio, ro), for the following values of the parameters: (a) ojx — — 1.85, Sl^ = 0.75, f = — ^o;(b) As 
before, but for v = 0; (c) The same, but for v = +vq, (d) lux = —0.93, $7° = — 2, v = 0.96. 



For experiments like SNAP jjj] the relevant redshift range (a window of small and 
moderate redshifts up to z = 1.7) is far away from the asymptotic regime. Here the 
rich panoply of behaviors exhibited by the effective EOS of the AXCDM model is very 
evident in e.g. Figs. || and ||. Some of these traits should in principle be observable in these 
experiments and should be the most conspicuous ones. However, as repeatedly emphasized, 
the effective EOS behavior of the AXCDM model in that redshift range can approximate 
to a great extent that of the ACDM model - either because it is very close to uj e = —1 
or because there is a gentle transition (from phantom-like behavior into quintessence one) 
while keeping an average value oj e = — 1. 

7.5 Evolution of the ratio r = po/Pm and the cosmological coincidence problem 



In Fig.[T(J we display the ratio r, Eq. fl4.1| ), in units of the current ratio (ro — 7/3) as a 
function of the cosmological redshift variable, for various values of the parameters. All the 



plots in Figs. 1C and 12 exhibit a maximum below 10. However, let us clarify that there are 



regions of parameter space where the maxima can be much higher, but in any case finite. 
Having a model where the ratio r = po/ Pm is bounded for the whole Universe lifetime 
provides a natural explanation for the cosmological coincidence problem. One may view 
the property ro ~ 1 as a simple corroboration of the fact that we live at a cosmic time near 
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Figure 11: Comparison of the maximum of the ratio r = po/Pm f° r ^x 
in Fig. fLOa) and different values of v. 
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the transition redshift Zt where the Universe changed from decelerated expansion into an 
accelerated one. But this observation leaves the coincidence problem in the same place. For 
instance, in the standard ACDM model, once the redshift z% is trespassed by the evolution 
the ratio r(z) starts to grow without end. So it is a kind of puzzling "coincidence" that we 
just happen to live so near the transition point Zt- This assertion should not be understood 
on anthropic grounds. When we say that it is a coincidence to "live near the transition 
point zt" we mean that the probability to find a lapse of the Universe's lifetime where 
r(z) ~ 1 should be essentially zero for a Universe (like the standard ACDM) where this 
ratio is larger than any given number if we await enough time. In the AXCDM model, 
instead, the situation is quite different. If the Universe would still be populated by some 
form of intelligent life in the remote future, say 5 — 10 billion years from now, and these 
beings would be curious enough to perform a new measurement of r(z) at their time, they 
might learn that r(z) is no longer ?"o, as it is now, but maybe something like 2r$ or 5ro, or 
even back to tq at some later time (where the Universe would in this case be decelerating 
rather than accelerating). There would be no correlation between the value r(z) ~ 1 and 
the state of acceleration of the Universe, and they would not be compelled to conclude that 
they are living in a very peculiar spell of time. We can confirm this picture by looking at 



some numerical examples in Fig. |10j and |llj. In particular, in Fig. |11| we demonstrate the 
dependence of the maximum on the value of v for fixed values of the other parameters. 

For a more physical grasp of the situation, we have also solved numerically the AXCDM 
model in terms of the cosmic time variable t (cf. Fig. 12). In this case it is not possible 
to provide analytic expressions such as (|3.8|) or ( |5.7|) in terms of the z and £ variables 
respectively. We limit ourselves to provide the numerical result for the expansion rate 
H = H(t) (Fig.|l2|a) and the ratio r = r(t) (Fig. 12 b) in units of the present ratio tq for the 



cases (a),(b) and (c) of Fig. 10. We can see that in terms of the cosmic time the distribution 
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of the relevant events in the future (such as the reaching of the maximum and the arrival 
at the turning point) is not as sharp and abrupt as it was in terms of the redshift variable. 
Indeed, one can see that in terms of the cosmic time the maximum takes place far away in 
the future, typically at around 2 Hubble times (Hq ) since the beginning of the expansion, 
i.e. roughly one Hubble time after our time, where Hq ~ \AGyr for ho = 0.7. The 
stopping point can typically lie at around t s ~ (2.5 — 3) Hq . To avoid cluttering, scenario 
(d) of Fig. |l0| is not included in Fig. [l^, but we have checked that it amounts to a lower 
maximum (r < 4), and the stopping point lies at t s > 3.5 Hq . 

From our detailed analysis of the ratio r = pd/ Pm we may assert that the AXCDM 
model provides a possible solution (and in any case highly alleviates) the cosmological 
coincidence problem. At the very least it offers an alternative explanation to some other 



interesting proposals discussed in the literature [41, [4^, [4^, [55|]. For example, in phantom 
dominated universes one may compute the fraction of the lifetime of the Universe where the 
ratio (|4.3| ) stays within given bounds before the doomsday [43|. Recall that these universes 



are threatened and eventually killed by a final Big Rip, so that the idea is to compute the 
time interval between the start of acceleration and the time when the Universe is dismantled 
to its ultimate fundamental components, and compare it with the total lifetime. If the 
fraction is not very small, one may feel satisfied with this explanation. However, we can 
see that a purely phantom Universe is just a very particular case of AXCDM Universe, 
namely it corresponds to Ao = v = and ax < 0, Q x > 0. In a general AXCDM Universe 
the range of possibilities is much larger, e.g. there can be a turning point rather than a 
Big Rip, the ratio ( [4.3|) can stay bounded during the entire span of the Universe lifetime, 
and finally the matter is not destroyed. 

Recall from the analytical discussion in Sections ^| and || that we have three critical 
points in the future: z s (stopping/turning point), z$ (point where the total DE density 
vanishes and uj e becomes singular), and z max (point in the future where the ratio r(z) 
attains its maximum). They are ordered as follows: 

z s < z < z max . (7.20) 



In correspondence we have the following values of the ratio (4.1): 



r(z s ) = -1< r(zo) = < r(z max ) < 10 . (7.21) 

At the same time we may compare the evolution undergone by r(z) from the nucleosynthesis 
time through the present time till the maximum. We have exemplified typical situations 
in which 

\r N \~0.1<r ~l<r(z max )<10. (7.22) 

Namely, the ratio of the DE density to the matter density undergoes a slow evolution 
within a few orders of magnitude. In this context it does not look bizarre that its present 
value is of order one. It should, however, be clear that the condition po/Pm ~ 1 is 
not guaranteed for any given point of the AXCDM parameter space. In this paper we 
have not addressed the computation of the probability distribution for the two energy 
levels pd and p m to be comparable at a particular time of the cosmic evolution. In the 
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AXCDM model discussed here the complete solution, including the dynamics of the energy 
density levels at a particular time, depends deterministically on the set of values of the 
free parameters in Eq. ( |7,3|) . If we would discretize the parametric space, i.e. dividing 
it into small 3D elements, and then calculate the dynamics for every volume element, it 
would be possible to numerically construct the joint probability distribution for the matter 
and (composite) dark energy density levels at some fixed time. Since the volume of the 
parametric space is infinite, one would need to introduce additional constraints in order 
to select a finite volume of parametric space and make this calculation useful. Although 
an explicit construction of this distribution could indeed be performed after defining these 
constraints, our presentation actually provided an alternative procedure which should serve 
equally well our purposes; namely, we explicitly computed that part of the total parameter 
space which results in the wanted dynamics (that is to say, the ratio r = po/Pm being 
bounded from above with the value of maximal r within an order of magnitude from the 
present one) and showing that it is non-negligible. This is patent in Fig. |2[ where the 
"physical region" (viz. the one satisfying the aforementioned requirements) is exhibited. 
What we have found is that the dimensionality of the " physical region" is the same as the 
dimensionality of the entire parametric space itself, in other words: the "physical region" is 
of non-zero measure in the space (fi^> wj, v). Not only so, the fractional volume (referred 
to the volume of the region where the parameters stay within their "natural" values, see 
Fig.g) is sizeable, as shown by the various projections onto the different planes displayed 
in Fig. y-Q. This proves that no additional fine-tunings are necessary to realize this 
scenario which is, in that sense, natural within the AXCDM models. This fact, combined 
with the nice features of the effective EOS of dark energy (described in Section |7.2| ) at 
the redshifts amenable to distant supernovae observations (and the compatibility with the 
nucleosynthesis bounds), further reinforces the AXCDM models as models with acceptable 
behavior of dark energy which at the same time provide an acceptable explanation (or 
substantial alleviation) of the coincidence problem. 

7.6 Transition from deceleration to acceleration 

A precise measurement of the transition redshift zt from deceleration into acceleration 
could also help to restrict the parameters of a given model. The value of zt is obtained 
theoretically by considering the change of sign of the deceleration parameter from q{z) > 
to q(z) < 0, where 

q(z) = ~ = a = - \{1 + 3w m ) n m (z) + (1 + Zuo e ) n D {z)] . (7.23) 

The cosmological functions (Id(z) and Q m (z) are defined in ( |3.14| ). Since the the point z% 
is reached during the matter epoch, where oj m = 0, tu e < 0, we can define zt as the point 
that satisfies q(zt) = 0, hence 

S£M = , \ . (7.24) 

^m(zt) 3|w e |-l 

In the standard ACDM model the corresponding relation reads Vto{zt) = Q m (zt)/2, which 
is recovered from ( 7.24| ) after setting uj e = — 1 and identifying tip with O^. Equivalently, 
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Figure 12: (a) The expansion rate of the AXCDM model in terms of the cosmic time, H = H(t), 
for the same parameter choices as in Fig. |lO|a,b,c; (b) The ratio r — po/pm, in units of its present 
value ro, as a function of the cosmic time for the same scenarios as in (a). In both cases the cosmic 
time is measured in Hubble time units H^ 1 . The present time lies at to = 13.7 Gyr, i.e. ~ 0.98 -ffo" 1 



the deceleration parameter can be computed from the expansion rate H(z) as follows 

1 ,. ,1 dH 2 (z) 



< {Z) = - 1+ 2 {1 + Z) HH Z ) dz 



(7.25) 



and zt satisfies q(zt) = 0. For the standard ACDM model and using the favorite values 



^M = O-^ an d ^a = 0-7 °^ the cosmological parameters in the flat case, one obtains 



z t 



-1 + 



f^H 67 - 



(7.26) 



For the AXCDM model the expansion rate H = H{z) is more complicated and is of 
course sensitive to the basic set of parameters ( [T.3| ) through the formula given in the first 
expression in Eq. (|3.8|). However, due to the relatively wide experimental margin obtained 
for zt (cf. Riess et al. in M) 



cxp 



0.46 ±0.13. 



(7.27) 



it is not possible at the moment to firmly discriminate which value of v is preferred. Notice 
that the range ( [7.27 ) excludes the standard model value ( 7.26 ) by a bit more than one 
standard deviation, but this is not very significant for the moment. We expect that this 



situation will substantially ameliorate in the future. In the meanwhile from Fig. 13 we test 
the sensitivity of the transition redshift zt to the value of v in the AXCDM. It seems to 
favor positive values of v. It is interesting to note that the zt sensitivity to v is opposite 



in the AXCDM as compared to the pure running A model of Ref. [23|. In the latter the 
values v < are favored (see particularly Fig. 4 of the second reference in |23|] ) . This is 
an example of how the observables can discriminate between two different realizations of 



a cosmological RG model. Both models, the AXCDM and the running A model of [23] 
are based on a similar RG framework. However, in the latter matter is not self-conserved 
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Figure 13: Comparison of the transition rcdshift z t (deceleration-acceleration) for lox 
0.75) and different values of v. 
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(there is exchange of energy between matter and the running A) whereas in the present 
model matter is self-conserved and the running A can exchange energy only with the 
cosmon. Obviously, it will be necessary a comparison of several high precision observables 
(a task presumably feasible thanks to the next generation of observational experiments) 
before we can discriminate between the two kinds of RG cosmologies. Another observable 
that should help to discriminate between models is the EOS itself and the existence of a 
transition point z* where the model might effectively change from quintessence to phantom- 



like behavior or vice versa. For example, the RG model of [23] may have transitions from 



quintessence into phantom, as was shown in Ref. [37], whereas the AXCDM model can only 
have transitions of the opposite kind, or no transition whatsoever if stopping is required 
(cf. Section ||). Let us also add that the RG model of [23] does not have any clue, 
in contradistinction to the present AXCDM model, as to finding an explanation for the 



cosmological coincidence problem because the ratio (4.1) - which in that model just reads 
Pa{z)/pm{z) - is unbounded. 

7.7 Phase trajectories of the cosmological autonomous system 



Finally, in Fig. |l4|and 15 we show some phase trajectories of the autonomous system ( |5.2[ ) in 
two different planes and for particular scenarios covering the two types of situations (node 
and a saddle point). Due to the stopping conditions discussed in Section[5], in both of 
these cases the trajectories become eventually truncated and therefore either did not have 
a chance to reach the convergence node, or alternatively fail to escape away to infinity, 
respectively. In Figs. |l4|a,b we display several phase trajectories in the (Qa,^x) plane 
covering these possibilities. In Fig.|l4|a we have an example that fits with scenario VII of 
Section ||. The trajectories that turn left and get stopped are those fulfilling the condition 



(5.23), whereas those turning right do not fulfill it and hence do not stop. In contrast, in 
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Figure 14: Phase trajectories of the autonomous system (5.2) in the (p,\,Clx) plane. Dashed lines 
show the parts of the curves corresponding to positive redshifts (i.e. our past), full lines indicate the 
parts between the present moment and the stopping (if there is stopping) and dotted lines denote 
the inaccessible part of the trajectory after the stopping, (a) Example of scenario VII in Section^ 



(saddle point case) corresponding to u>x = —1.85, v = —vq. The eigenvector directions (5.5) are 
also marked. The various trajectories correspond to different choices of f2 A indicated in the figure; 
(b) Example of scenario IV with a node determined by {y, fi^)) m ^ ne n g ure ^a = — 2, v = 0.96. 
All trajectories are attracted towards the node ( 5.15| ), but they get stopped on the way to it. The 
point £ = corresponds to the present. The cosmic evolution is in the £ > direction, see Eq. (5.1). 
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Figure 15: Phase trajectories of the autonomous system (5.2) in the (0_D,fi m ) plane. Cases (a) 
and (b) are respectively as in Figs. pT|a,b but here we plot the corresponding trajectories in the 
phase plane for the total dark energy density versus matter density, i.e. (Qo,^m)- The meaning 
of the lines is as in Fig. [UL 



Fig. 14b we are dealing with scenario IV, where in the absence of stopping the trajectories 
would be attracted to the node ( 5.15|) for all values of u>x in the quintessence range (u>x ^ 
—1). In this particular example, which illustrates a situation where Q\ and v are both 
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rather large, the fixed point is well out of the figure. Even so the trajectories show already a 
clear trend to get focused. The choice of the parameters does satisfy the stopping condition 



( 5.16 ) and so the trajectories get stopped well before reaching the node. It can be clearly 



seen in the figure because the "stopped part" of the trajectories is partially shown. 



In Figs. 15a, b we represent respectively the same kind of scenarios as in Fig. 14a,b 
except that we perform the representation in the plane of total DE density versus mat- 
ter density. In Fig.|i~5|a the trajectories flp satisfying the condition ( |5.23 ) are seen to get 



stopped when crossing the straight line Qd = —&-AI corresponding to r = — 1. The tra- 
jectories that do not fulfill ( |5.23 ) escape away into unbounded values of &d- Similarly, 



in Fig.|lqb we have a situation with the presence of a "would be node" which focus all 



trajectories, but all of them get stopped because the condition fl5.16|) is fulfilled. 

8. Discussion and conclusions 

In this paper we have explored the possibilities offered to cosmology by a non-trivial com- 
posite model of the dark energy (AXCDM). By composite model we mean that it contains 
at least two entities of entirely different nature, i.e. not just a replica of a given one (e.g. 
two or more scalar fields). By non-trivial we mean that these two components generally in- 
teract with one another. From our point of view it is perfectly conceivable that the observed 
DE is made out of a mixture of various components (perhaps some scalar fields) stemming 
from different sources of a fundamental theory, e.g. string theory. However, in the most 
general case we have to contemplate the possibility that one of the ingredients of the DE 
could be a cosmological term, A, with its own dynamics and interactions with the rest. In 
our model the DE components other than A have been represented by X, the "cosmon". 
We should stress that the cosmon need not be a fundamental field, it could be an effective 
representation of dynamical fields of various sorts, or even the effective behavior of higher 
order curvature terms in the effective action. In this paper we have presented a composite 
DE model of this kind, which we have narrowed down to its simplest non-trivial form. The 
dynamics of the variable cosmological term, A, is linked to the renormalization group (RG) 
running in QFT, and the effective entity X interacts with A, but not with matter, such 
that both the total DE and matter-radiation densities are separately conserved. In fact, 
the cosmon (whose ultimate nature is not necessary to specify in our discussion) follows an 
evolution which is fixed by the RG dynamics of A itself through the requirement of con- 
servation of the total DE. The particular modeling of the A evolution that we have used 
follows the law 5 A ~ H 2 , which is well-motivated in the RG framework of Ref. |2^, 23, 45]. 



The modeling could of course be different in other implementations of the RG in Cosmo- 
logy, or even without explicit reference to it. In this sense the class of AXCDM models 
is surely more general. However, the non-trivial implementation presented here furnishes 
a "proof of existence" of promising scenarios of this kind, and suggests that the idea of 
having a "scaling A" may constitute a fruitful interface between QFT and Cosmology [£0[ . 
At the very least it gives a hint on the immense capability of the AXCDM cosmologies to 
describe the next generation of high precision data, even without making explicit use of 
simple scalar field models of the DE with particular forms of their potentials. 
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In the AXCDM context the A term has a very different status as compared to the 
ACDM model. While the main observable quantity here is the total DE density po = 
PA + Px, the AXCDM model places A in a fully dynamical context together with the 
cosmon. Within this approach the puzzling question, formulated within the standard 
ACDM model, of why A takes a particular value at present, and why this value is so 
close to the matter energy density does no longer apply. The total DE is an evolving 
quantity through the entire history of the Universe, and under appropriate conditions it 
traces closely the matter density. This is so because in this model we can find a large 
region of the parameter space where the expansion of the Universe stops at some point in 
the future. In contradistinction to the standard ACDM model, here the ratio pd/ Pm is 
not divergent with the evolution, and it reaches a maximum at some point in our future 
before the Universe stops and reverses its motion. We may view the present state of our 
Universe in this model as just being in transit to reach that maximum and the subsequent 
stopping/turning point. At present that ratio is of order 1, but at the maximum it can be, 
say ten times bigger. 

Obviously, this framework provides some clue for understanding the cosmological co- 
incidence problem. From our point of view this problem is not so much the issue of why 
Pd/ Pm is of order one at present, but rather of why po/Pm — ► oo in the asymptotic future 
regime (as within the standard ACDM model). If that ratio is bounded, as it is the case in 
the model under consideration, then the Universe can be at any value below this maximum, 
and this feature certainly improves the status of the purported cosmological coincidence. 
Within the AXCDM model the fact that pd/ Pm ~ 1 can be looked upon as a reflex of 
the present state of accelerated expansion of the Universe, which started somewhere be- 
tween z = 0.3 — 0.7, cf. Eq. ( 7.27| ). This interpretation is not possible in the standard 



ACDM model because once the situation pn/Pm ~ 1 is reached the ratio further increases 
without end, and so it is hard to understand why at present this ratio comes out to be of 
order one. It can only be if we admit that we live in a very special moment, namely one 
rather close to the start of the acceleration. Similar problems are met by the quintessence 
and phantom models in the literature. In the former case we cannot understand why the 
(slowly ever-decreasing) quintessence field matches up the (fast decreasing) matter density 
precisely right now in the context of an eternal Universe, whereas in the latter we cannot 
figure out why the (slowly ever- increasing) phantom field just catches up with the matter 
density squarely around our time. 

Quite in contrast, in the AXCDM model, with non-vanishing v, the ratio pd/ Pm re- 
mained relatively small in the remote past, being roughly proportional to \v (l+wx)j. (We 
can easily provide values of the parameters for which po/Pm ~ 10% at the nucleosynthesis 
era.) Later on, it started to be of order 1 around the transition epoch (z = 0.3 — 0.7) where 
the Universe changed from a decelerated expansion into an accelerated one. From the early 
times up to the present day the ratio increased slowly, and in the future it will continue 
increasing, but it will never surpass a finite number, e.g. of order 10, until the Universe 
will stop and reverse its motion (a few Hubble times after our time). The upshot is that 
from the time of primordial nucleosynthesis till the remote point in the future where the 
ratio pd/ Pm becomes maximum, it changed only mildly, say from 0.1 to 10, being currently 
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around 1, in contrast to the fantastic change predicted in the ACDM case (viz. from 
to oo) during its infinite lifetime. Therefore, once the AXCDM Universe starts to acceler- 
ate, the probability to live at a time where r is of order 1 is essentially 1. Obviously, this 
model can significantly ameliorate the situation with the cosmological coincidence problem, 
and it may actually serve as a prototype example of how to reach an explanation for this 
"coincidence" through the addition of some X to ACDM. 

Furthermore, we have shown that the effective EOS behavior of the AXCDM model 
can mimic to a high extent that of the standard ACDM model (namely by providing an 
average value of the effective EOS parameter uj e near —1) in essentially the entire redshift 
interval where the SNAP experiment M is going to search for distant supernovae in the near 
future. Moreover, the possibility of interaction between the cosmon and the cosmological 
term (which is gauged by the parameter v) allows to modulate the effective EOS of the 
AXCDM model by giving it a tilt into the quintessence or phantom regime depending on 
the sign and size of v and the value of the barotropic index of the cosmon, ux ■ With no 
interaction (y = 0) between the two components of the DE, this would be impossible. This 
shows that even though Q\ and Qx may not be individually observable, the interaction 
between A and X can be detected through the effective EOS features that it gives rise to. 

We wish to stand out another important signature of the model that involves compar- 
ing the two main ranges of redshift observations. Remarkably, while the behavior of the 
effective EOS of the AXCDM model at small and intermediate redshifts (relevant to the 
supernovae analysis) can be very close to the standard ACDM model, it turns out that at 
the very high redshift range (relevant for CMB considerations) the effective EOS of the DE 
is similar to that of matter-radiation. In this range the DE physical effects of the AXCDM 
model could be detected through a "renormalization" correction of the cosmological mass 
parameter f2^ (which could be up to order 10%) as compared with the value fitted from 
the supernovae data. It turns out that most of the parametrizations used in the literature 
to (simultaneously) fit the effective EOS of the DE at low and high z could be inadequate 
to describe the AXCDM model and other promising RG cosmological models whose EOS 



behavior exhibits this kind of features [37]. 



We may wonder if an effective description of the composite dark energy used in the 
AXCDM model can be performed in terms of two or more interacting scalar fields. The 
problem essentially boils down to finding an appropriate potential for the scalar fields. 
Should the form of this potential be a very complicated or a contrived one, the effective 
description in terms of scalar fields would not be very appealing. A simple example shows 
that the scalar field description is not always attractive. Take the special case of the 
AXCDM model in which v = 0, i.e. the dark energy components do not mutually interact; 
in particular, we could think of one component being a negative A and the other being 
e.g. a quintessence component X. In this case an effective description in terms of two 
non-interacting scalar fields would not be especially useful since one of the fields would 
have to be non-dynamical. If on the other hand they would be taken as interactive, this 
would lead to an unnecessarily more complicated scalar field description of the DE as 
compared to the original AXCDM one, and therefore it could not be considered as a better 
option. Furthermore, even if an effective scalar field model is found which reproduces the 
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global expansion of the universe in the AXCDM model, local properties such as growth of 
perturbations might be different in these models. As we stressed before, the cosmon entity 
need not be a fundamental field. In this sense the AXCDM model can be the effective 
behavior of many possible constructions, which in all cases have in common the existence 
of a conserved total DE density. This warns us that we should be careful to avoid extracting 
too precipitated conclusions on the ultimate nature (e.g. scalar field) of the DE. 

Clearly, the AXCDM model can encompass very general behaviors of the DE. There- 
fore, the fact that the standard ACDM with a strictly constant A seems to be in fairly good 
agreement with the last round of high precision WMAP experiments ||] - with the possi- 
bility that future SNAP data could independently corroborate this result - in no way can 
be considered a definitive confirmation of the ACDM. As shown in this paper, the AXCDM 
model could do in principle a very similar job, with the bonus that it can offer a natural 
explanation for the cosmological coincidence problem. It remains to study its implications 
for the structure formation and CMB, but this is of course a task beyond the scope of 
the present paper. For the moment we have shown that the model can be perfectly tested 
by performing measurements at all accessible ranges of cosmologically significant redshifts 
(both low and high), and therefore we expect that the next generation of high precision 
experiments [j|], ||, ||] will allow to better pin down the features of the effective EOS of the 
dark energy such that the AXCDM parameters can be better determined. 

We conclude with the following observation. The "addition" of a dynamical cosmon 
X to the standard ACDM model can be viewed (or not) as a natural feature in DE model 
building. But in practice a cosmon-like entity (e.g. in quintessence regime) has thrived 
alone through the literature quite successfully (as a full substitute for A) for quite some 
time, following perhaps the courageous - and certainly longer - tradition of A itself. On 
the other hand, one may ponder whether it is natural or not to allow a dynamical A term. 
However, the running of A, or of any other parameter, may be viewed as a natural property 
of Quantum Field Theory applied to Cosmology (2(]] . At the end of the day it seems that 
other cosmological problems could be significantly smoothed out if the DE is permitted to 
be neither pure A nor pure X, but just some more or less democratic combination of both, 
including some peaceful interaction between them. 
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